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I .  INTRODUCTION 


A  rash  of  railway  accidents  involving  tank  cars  carrying  flammable 
materials  such  as  propane  and  vinyl  chloride  prompted  an  investigation 
by  the  Federal  Railroad  Administration  (FRA)  of  the  Department  of  Trans¬ 
portation.  Because  of  the  Ballistic  Research  Laboratory's  (BRL)  knowl¬ 
edge  in  fire  technology  and  skills  in  high-temperature  measurements,  the 
FRA  sought  out  the  BRL  for  assistance.  A  series  of  field  experiments  on 
one-fifth-scale  model  tank  cars  culminated  in  two  full-scale  pool  fire 
tests  on  rail  tank  cars  filled  with  propane. 

For  the  first  test,  a  bare  car  was  placed  in  a  pit,  instrumented 
and  filled  with  approximately  113.5  kl  (30,000  gallons)  of  propane.  The 
car  was  then  completely  engulfed  in  flames  (JP-4  jet  fuel  was  used  to 
provide  the  source  of  the  fire).  After  24.5  minutes  of  fire  exposure, 
the  tank  violently  ruptured,  spewing  fire  and  debris  considerable  dis¬ 
tances.  Armed  with  the  data  from  the  test,  post-test  analyses  of  the 
fragments,  metallurgical  analyses,  etc.,  the  mechanism  of  failure  was 
understood  and  documented. I 

Propane  and  vinyl  chloride  are  shipped  under  pressure  in  the  liquid 
form.  Tn  a  fire  environment,  heat  is  transferred  to  the  interior  of  the 
car.  The  temperature  of  the  contents  is  raised  and  the  pressure  in¬ 
creases.  Above  a  specified  pressure,  the  relief  valve  opens  and  vents, 
but  the  pressure  will  continue  to  climb  if  the  total  heat  to  the  interi¬ 
or  is  great  enough.  In  addition,  the  steel  tank  shell  is  heated,  par¬ 
ticularly  that  portion  which  is  not  in  contact  with  the  liquid.  As  the 
metal  heats,  its  strength  begins  to  degrade,  especially  where  metal  tem¬ 
peratures  are  in  excess  of  425°C  (800°F) .  If  exposed  to  the  fire  long 
enough,  stresses  induced  by  the  interior  pressure  exceed  the  strength  of 
the  temperature-weakened  steel.  A  sudden  release  of  the  pressurized 
contents  then  occurs,  vaporizes,  mixes  with  air,  and  an  explosion  occurs. 
This  fuel/air  explosion  can  "rocket"  portions  of  a  tank  car  for  consid¬ 
erable  distance,  create  a  tremendous  fireball  and  a  damaging  blast  wave. 

The  second  full-scale  test  was  on  an  insulated  car,  and  it  was  dem¬ 
onstrated  that  by  protecting  the  metal  skin  from  reaching  elevated  tempera¬ 
tures,  the  car  would  retain  its  integrity  while  releasing  its  contents 
through  the  relief  valve  until  the  car  was  essentially  empty.  Also,  the 
rupture,  if  indeed  it  did  rupture,  would  be  significantly  less  violent." 


Charles  Anderson ,  William  Townsend ,  John  A.  Zook,  Gregory  Cowgill,"The 
Effects  of  a  Fire  Environment  on  a  Rail  Tank  Car  Filled  With  LPG, ", BRL  R 
1935,  USA  Ballistic  Research  Laboratories,  Aberdeen  Proving  Ground ,  MD , 
September  1976.  (AD  B015605L) 
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Charles  Anderson,  William  Townsend,  John  Zook,  Gregory  Cowgill,"The 
Comparison  of  Thermally  Coated  and  Uninsulated  Rail  Tank  Cars  Filled  With 
LPG  Subjected  to  a  Fire  Environment ,  "  Report  No.  FRA-OR8  D  75-32,  National 
Technical  Information  Service,  Springfield,  VA,  December  1974. 


However,  accident  investigations  showed  that  a  significant  percent¬ 
age  of  ruptures,  where  there  were  reliable  estimates  for  the  times,  oc¬ 
curred  in  less  than  24.5  minutes.  Yet,  full-scale  testing  of  the  base  car 
was  a  full  engulfment  pool  fire  -  a  "worst  case"  test.  Hence,  another 
mechanism  was  identified  as  contributing  to  tank  car  failures  -  the 
"torch"  mechanism. 

A  torch  can  result  from  ignited  liquid  and/or  vapor  rushing  out  of 
a  hole  or  tear  in  the  shell.  This  could  be  caused  by  a  coupler  impact 
in  an  accident,  or  possible  effluent  from  a  relief  valve  on  an  overturned 
car  impinging  on  an  adjacent  car.  Remembering  that  the  contents  are 
already  under  pressure,  and  that  any  additional  heat  from  a  fire  would 
raise  the  internal  pressure,  the  ignited,  high  velocity  stream  acts  as 
as  large  blowtorch.  High  heat  fluxes  from  the  torch  can  heat  the  im¬ 
pinged  steel  to  dangerously  high  temperatures  in  only  several  minutes. 

To  investigate  this  phenomenon,  the  BRL  designed,  constructed,  and 
operates  a  torch  simulator  facility  at  the  DOT  Transportation  Test  Center, 
Pueblo,  Colorado. There,  the  BRL  tests  the  thermal  response  of  steel 
plates  with  various  types  of  thermal  insulating  systems  for  possible  ap¬ 
plications  to  rail  tank  cars. 4  The  torch  is  formed  by  the  blowdown  of 
pressurized  propane  vapor  and  liquid  from  a  1.9  kl  (500  gallon)  tank. 
Figure  1.  The  vapor  and  liquid  are  combined  in  a  manifold  prior  to  ex¬ 
pansion  through  a  0.95  cm  (3/8-in)  orifice.  The  resulting  high-velocity 
jet  is  initially  ignited  and  the  resulting  torch  (a  turbulent  diffusion 
flame)  impinges  on  the  test  plate. 

Analyses  of  the  torch  velocity,  its  downstream  lateral  spreading, 
and  the  heat  flux  to  the  plate  are  useful  for  several  reasons.  Initially, 
BRL  needed  guidelines  to  help  design  the  tests.  Then,  after  demonstrating 
that  a  torch-like  test  is  a  reliable  way  to  characterize  the  performance 
of  thermal  coatings,  BRL  wanted  to  understand  the  effects  of  various  pa¬ 
rameters  and  learn  what  scaling  parameters  characterize  the  torch.  For 
small-scale  tests,  an  analysis  is  needed  to  specify  the  torch  parameters 
that  adequately  simulate  large-scale  heat  fluxes.  Southwest  Research 
Institute,  therefore,  was  requested  by  BRL  to  analyze  the  fluid  dynamics 
and  thermodynamics  of  a  high-speed  propane  torch  impinging  upon  a  plate, 
to  include  propane  jet  formation  and  combustion. 


William  Townsend  and  Richard  Markland,  " Preparation  of  the  BRL  Tank  Car 
Torch  Facility  at  the  DOT  Transportation  Test  Center,  Pueblo,  Colorado, " 
BRL  unp  ublished  report,  USA  Ballistic  Research  Laboratories ,  Aberdeen 
Proving  Ground,  MD,  September  1975. 

4 

Charles  Anderson,  William  Townsend,  Richard  Markland,  John  Zook,  "Compari¬ 
son  of  Various  Thermal  Systems  for  the  Protection  of  Rail  Tank  Cars  Tested 
at  the  FRA/BRL  Torching  Facility , "BRL  unpublished  report,  USA  Ballistic 
Research  Laboratories ,  Aberdeen  Proving  Ground,  MD,  December  1975. 


Figure  1.  Schematic  Representation  of  Two-Phase  Torch  Test 


II.  JET  ANALYSIS5 

At  the  exit  plane  of  the  propane  tank  orifice,  a  two-phase  propane  jet 
emerges  at  supersonic  velocities  as  the  high  pressure  fuel  expands  to  atmos¬ 
pheric  pressure.  The  jet  then  spreads  laterally  by  the  turbulent  entrainment 
of  air.  Ignition  can  occur  where  sufficient  air  has  been  entrained  to  support 
combustion.  At  this  point,  a  new  jet  is  formed  by  the  products  of  combustion 
and  continues  to  spread  laterally  as  more  air  is  entrained  through  the  turbulent 
mixing  process.  As  the  jet  nears  the  target  plate,  the  flow  is  turned,  which 
causes  a  decrease  in  the  centerline  velocity  of  the  jet,  and  a  wall  jet  is  formed. 

A  rigorous  two-dimensional  treatment  of  this  prgblem  may  be  possible  by 
using  techniques  such  as  those  reported  by  Spalding.  However,  for  speed  and 
efficiency, a  simple  analysis  has  been  developed  by  combining  a  number  of  one¬ 
dimensional  solutions  available  in  the  literature.  This  approach  is  relatively 
straightforward,  physically  sound  and  sufficiently  accurate  to  explain  observed 
experimental  phenomena. 

A.  Reservoir  Expansion 

The  propane  torch  analysis  is  initiated  by  establishing  the  state  of  the 
fuel  upstream  of  the  orifice.  Referring  to  Figure  1,  it  is  assumed  that  a 
mixture  of  saturated  vapor  and  saturated  liquid  occurs  in  the  manifold  upstream 
of  the  orifice.  The  pressure  is  monitored  at  this  point  during  the  test; 
however,  this  is  insufficient  to  establish  the  thermodynamic  state  of  the 
mixture.  Since  the  propane  is  a  mixture  of  both  vapor  and  liquid,  the  relative 
amounts  of  these  two  phases  must  be  established  in  order  to  determine  the 
thermodynamic  state  of  the  reservoir.  This  is  accomplished  by  defining  the 
"quality"  of  the  mixture  as  the  ratio  of  the  mass  of  vapor  to  the  total  mass  of 
the  mixture.  A  quality  of  one  represents  a  pure  vapor  flow,  while  a  quality 
of  zero  represents  a  pure  liquid  flow.  It  is  possible  to  control  the  quality 
of  the  mixture  by  a  suitable  adjustment  of  valves  A  and  B  shown  in  Figure  1. 

For  the  purpose  of  this  analysis,  the  reservoir  quality,  XR,  will  be  assumed 
to  be  known. 

The  thermodynamic  properties  of  propane  are  known  and  have  been  inserted 
into  this  analysis  in  tabular  form.  The  reservoir  temperature  is  found  by 
linear  interpolation  of  the  reservoir  pressure  since  the  temperature  is 
independent  of  quality.  The  remaining  state  variables,  such  as  entropy, 
enthalpy  and  specific  volume,  are  then  found  on  the  basis  of  the  assumed 
quality. _ 

C 

tJA  more  complete  description  of  the  modeling  is  given  in  "An  Analytical 
Model  of  a  Two-Phase  Propane  Combusting  Jet,"  James  L.  Rand,  Report  02- 
5045-001,  Southwest  Research  Institute,  San  Antonio,  TX,  November  1980. 

tj Spalding ,  D.B.,  "A  Simple  Model  on  the  Rate  of  Turbulent  Combustion ," 

Turbulent  Combustion, Progress  in  Astronautics  and  Aeronautics,  Vol.  58, 

AIM,  1978,  pp. 105-116. 


16 


Cl  -  a) 


SR  =  (i  -  X)  SFR  ♦  (X)  SGR 
hR=  a  -  x)  hFR+  (X)  hGR  ci.b) 

vR  =  (i  -  x)  vFR  ♦  (X)  vGR  (l.c) 


rhe  thermodynamic  properties  on  the  right-hand  side  of  these  equations  are 
determined  from  the  tabular  data  in  the  same  manner  as  the  temperature. 

After  determining  the  reservoir  conditions,  the  fuel  is  then  expanded 
from  the  measured  high  pressure^through  Jh0.95-cm  orifice  to  ambient 
pressure.  In  a  previous  report  the  fuel  was  considered  to  be  an  ideal 
gas  which  permitted  the  jet  throat  and  exit  velocities  as  well  as  state 
variables  to  be  determined  analytically.  However,  due  to  the  saturated 
nature  of  the  reservoir,  a  numerical  technique  was  developed  which  would 
permit  the  expansion  of  the  propane  along  an  isentrope  through  the  orifice 
to  the  final  jet  exit  conditions. 

The  procedure  used  to  obtain  the  expanded  fuel  conditions  is  based  on 
holding  the  entropy  determined  in  equation  1  invariant  as  the  pressure  is 
reduced.  If  the  flow  is  assumed  to  be  choked  at  the  orifice,  then 
conditions  at  the  throat  will  result  in  the  maximum  flow  rate.  Since  the 
effective  area  of  the  orifice  is  relatively  insensitive  to  small  changes 
in  pressure,  the  maximum  flow  rate  will  coincide  with  the  conditions  for 
maximum  mass  flux,  PTUT*  The  density  at  the  throat,  pT,  is  found  by 
incrementally  reducingthe  pressure  and  computing  the  quality  required  to 
maintain  the  same  entropy.  Having  determined  the  quality,  the  enthalpy,  hT, 
may  be  evaluated  and  the  flow  velocity  computed  from  the  energy  equation: 

hT  *  t  UT  -  hR  (2) 

This  process  is  repeated  until  the  pressure  for  a  maximum  value  for  the 
mass  flux  is  found.  Once  this  is  determined  all  other  thermodynamic 
properties  may  be  found  at  the  throat.  The  mass  flow  is  determined  from 
the  product  of  the  effective  area  of  the  throat  and  the  mass  flux. 

m  =  pTUTATCD  (3) 

The  coefficient,  Cp,  may  be  considered  to  be  a  discharge  coefficient  which 
accounts  for  the  vena  contracta  at  the  orifice.  The  value  of  this  coefficient 
may  be  adjusted  on  the  basis  of  experimental  observations;  however,  for  a 
ratio  of  orifice  diameter  to  inlet  pipe  diameter  less  than  0.2,  a  value  of 
0.5  appears  to  be  applicable  over  a  large  range  of  Reynolds  numbers. 


7 

Astleford ,  W.J.  and  Dodge ,  F.T. ,  "Response  of  Fuel  Targets  to  Munitions: 
Heat  Transfer  Analysis  of  Torches  Impinging  on  Plates , "  SWRI  Report  HO. 
02-3669 ,  June  1977. 
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The  exit  conditions  from  the  orifice  may  be  found  by  further  expansion 
of  the  propane  to  ambient  pressure.  If  the  exit  conditions  are  assumed  to 
be  one  atmosphere,  then  the  fuel  will  exit  at  a  temperature  of  231°K.  The 
quality  of  the  flow  may  be  found  from  equation  1  where  the  values  of  Sp^ 
and  are  known  from  the  table  of  thermodynamic  properties  at  one  atmos¬ 
phere:  Once  the  exit  quality  is  established,  all  other  static  properties 
of  the  flow  may  be  computed.  As  before,  the  velocity  is  computed  from  the 
change  in  enthalpy  of  the  flow  from  the  reservoir  enthalpy.  Typical  re¬ 
sults  which  show  the  influence  of  reservoir  conditions  on  the  exit  velocity 
are  shown  in  Figure  2. 


The  exit  area  is  computed  on  the  basis  of  the  mass  flow  rate  previously 
computed  at  the  throat  of  the  nozzle. 


P„IJ  A  C 
a  _  T  T  T  D 

e  p  U 
e  e 


(4) 


At  this  point,  all  flow  conditions  are  known  for  the  fuel  as  it  is  in¬ 
jected  into  ambient  air  considered  to  be  at  rest. 


B.  Mixing  Model 

A  wealth  of  literature  exists  on  the  expansion  and  turbulent  mixing 
of  a  free  jet.  The  work  of  Kleinstein  °  and  Witze^show  remarkably  good 
agreement  with  nonreacting  jLpt  mixing  experiments.  Portions  of  their  work 
were  used  in  a  prior  report  on  this  subject.  However,  in  order  to  sim¬ 
plify  the  analysis  here  and  maintain  the  one  dimensionality  of  this  tech¬ 
nique,  a  simplified  mixing  theory  has  been  developed  for  this  program. 

It  is  observed  that  as  the  fuel  exits  the  nozzle  no  air  has  yet 
been  mixed  with  the  fuel  so  that  combustion  cannot  occur.  The  flamma¬ 
bility  limits  in  air  of  propane  are  from  2.1%  to  10.1%  by  volume.  There¬ 
fore,  as  the  fuel  jet  expands,  the  mixture  will  reach  a  rich  flammable 
condition  first  and  then  approach  the  stoichiometric  ratio  further  down¬ 
stream.  In  addition,  the  flow  velocity  of  the  fuel-air  mixture  must  be 
slow  enough  to  sustain  combustion.  In  the  following  analysis,  the  fuel 
and  air  will  be  assumed  to  be  thoroughly  mixed,  which  is  a  departure 
from  those  assumptions  contained  in  References  8  and  9,  bo*  should  be 
sufficient  to  allow  an  estimate  of  the  products  of  combustion. 

It  is  realized  that  there  exists  conditions  in  which  there  is  an 
unburned  core  of  fuel-rich  propane,  with  burning  on  the  periphery.  As 
the  initial  conditions  go  to  more  and  more  liquid  (i.e.,  the  quality  goes 
toward  zero),  it  is  expected  that  this  fuel-rich  core  condition  will  be¬ 
come  more  prevalent  (and  indeed,  this  is  what  is  observed).  The  modeling 
assumption  is,  therefore,  not  as  well  met  as  when  the  quality  is  closer  to 
one. _ 

\leinstein,  G. ,  "Mixing  in  Turbulent  Axially  Symmetric  Free  Jets,"  J. 
Spacecraft  3  Rockets,  1,  July-August  OCA.  pp. 403-408. 

9  . 

Witze,  P. ," Center  line  Velocity  Decay  of  Compressible  Free  Jet,"  J, 
Spacecraft  &  Rockets,  12,  April  1974,  pp  417-418. 
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Figure  2.  Effect  of  Rese 


In  order  to  determine  the  relationship  between  the  velocities  and 
mass  ratios,  an  analytic  model  has  been  developed  on  the  basis  of  the 
conservation  of  mass  and  momentum  assuming  the  torch  does  not  change 
with  time.  The  steady  flow  situation  is  depicted  in  Figure  3(a).  Two 
regions  of  interest  will  be  considered  separately.  The  first  region  from 
the  orifice  exit  to  point  1  is  a  region  where  mixing  occurs  without  com¬ 
bustion.  The  second  region,  from  point  1  to  point  0,  is  where  instan¬ 
taneous  combustion  will  be  assumed  to  occur  without  the  addition  of  air. 

The  conditions  at  point  0  will  be  assumed  to  be  the  starting  conditions 
for  turbulent  mixing  of  a  free  jet  formed  by  the  combustion  products. 

A  control  volume  is  drawn  around  the  region  from  the  orifice  exit 
to  point  1  as  shown  in  Figure  3(b).  The  sides  of  this  control  volume  are 
parallel  to  the  centerline  and  sufficiently  far  away  that  the  velocity 
of  the  air  crossing  this  boundary  is  small.  It  should  be  noted  that  the 
velocity  of  the  pure  propane  fuel  crossing  the  left-hand  boundary  is  super¬ 
sonic.  Applying  the  principal  of  conservation  of  mass,  for  this  assumed 
steady  flow,  the  mass  flow  rate  of  the  mixture  is  given  by: 


0  A  U 
e  e  e 


(pAU) 


air 


=  PjAjU 


If  the  flow  is  steady,  then  the  ratio  of  the  mass  of  air  to  the  mass  of 
fuel  must  be  the  same  as  the  ratio  of  the  mass  rates  of  flow.  Therefore, 
the  mass  ratio  will  be  defined  as: 


m  = 


(oAU)  • 
_ air 

p  A  U 
e  e  e 


Thus  the  previous  equation  may  now  be  written  in  terms  of  the  known  mass 
flow  rate  of  fuel  and  the  mass  ratio  as: 


P1A1U1  =  (1  +  m)  PeAeUe  (5) 

Applying  the  principle  of  conservation  of  momentum  in  the  axial  direction: 

(pA)  +  [p  A  U  ]U  ♦  (pAU)  .  U  ..  =  (pA),  +  I  p. A.U  jll 

K  e  ye  e  el  e  'air  air  *  1  I  1  1  lj  1 

Since  the  pressure  at  all  points  on  the  boundary  of  this  control  volume 
is  equal  to  the  ambient  pressure,  the  external  force  terms  are  self  equil¬ 
ibrating.  The  axial  component  of  the  air  crossing  the  lateral  boundaries 
is  assumed  to  be  zero  by  the  judicious  location  of  the  boundary.  Therefore, 
this  equation  may  be  reduced  to  the  following: 


U 


1 


PeVM  « 

piVi)  e 


20 


(a)  Model  of  Jet  Starting  Condition 


(b)  Initial  Mixing 
Region 


(c)  Combustion  Region 


Schematic  of  Jet  Starting  Conditions 


Incorporating  equation  5  into  the  resulting  momentum  equation  yields  the 
velocity  of  the  mixture  in  terms  of  the  mass  ratio  and  the  propane  exit 
velocity. 

u  -  At 

U1  (1  +  m) 

Applying  the  same  reasoning  to  the  combustion  region  shown  in  Figure 
3(c),  conservation  of  mass  yields: 


p  A  U  = 
o  o  o 


PlAlUl 


(1  +  m)  p  A  U 
1  e  e  e 


Combining  this  with  conservation  of  momentum: 


U 


o 


U 

e 

(1  +  m) 


(6) 


The  area  over  which  the  combustion  products  are  moving  at  a  velocity  of 
Uq,  is  determined  from  the  mass  flow  rate.  Therefore, 


A 


o 


(7) 


In  order  to  compute  this  area,  it  is  first  necessary  to  evaluate  the 
density  of  tne  combustion  products.  It  should  be  noted  that  the  energy 
equation  has  not  yet  been  used. 


C.  Combustion  Process 


In  order  to  determine  the  thermodynamic  properties  of  the  combustion 
products  for  a  given  ratio  of  oxidizer  to  fuel,  a  program  developed  by 
NASA  Lewis  Research  Center^Owas  found  to  be  very  useful.  The  combustion 
of  propane  and  air  is  a  complex  chemical  equilibrium  process  which  re¬ 
sults  in  the  conversion  of  energy  into  a  sensible  form,  the  creation  of 
new  species  and  the  destruction  of  others.  The  program,  CEC76,  is  the 
latest  version  and  is  capable  of  computing  the  theoretical  thermodynamic 
properties  of  products  of  combustion  at  constant  pressure  and  energy. 


Gordon ,  S.  and  McBride ,  B.J. ,  "Computer  Program  for  Calculation  of 
Complex  Chemical  Equilibrium  Composition,  Pocket  Performance,  Incident 
and  Reflected  Shocks,  and  Chapman- Jouguet  Detonations , "  NASA  Lewis 
Research  Center  NASA  SP-273,  March  1976. 
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The  program  is  based  on  the  minimization  of  the  Gibbs  free  energy  rather 
than  the  use  of  equilibrium  constants  and  is  assumed  to  be  a  static 
process . 

In  order  to  relate  the  results  of  the  combustion  computations  to 
the  present  problem,  a  method  was  developed  to  account  for  the  various 
parameters  which  may  change  in  the  torch  test.  The  propane  fuel,  orig¬ 
inally  at  high  pressure,  will  be  injected  into  air  at  a  temperature 
corresponding  to  saturated  propane  vapor  at  one  atmosphere.  This  temp¬ 
erature  is  231°C  unless  the  flow  is  superheated  by  some  nonisentropic 
process  prior  to  injection.  The  air  is  assumed  to  be  at  rest  with  a 
pressure  and  temperature  of  one  atmosphere  and  298°K,  respectively.  The 
variables  which  may  influence  the  thermodynamic  properties  of  the  pro¬ 
ducts  of  combustion  include  the  quality  of  the  expanded  propane  jet  and 
the  energy  of  the  flow  which  may  be  represented  by  the  total  enthalpy  of 
the  jet.  In  addition  to  these  variables,  the  mass  ratio  of  air  to  pro¬ 
pane,  m,  plays  a  significant  role  in  the  efficiency  of  combustion,  the 
formation  of  species  and  the  resulting  flame  temperature  as  well  as 
other  thermodynamic  properties. 

The  NASA  program,  CEC76,  is  very  general  and  far  more  inclusive 
than  that  needed  for  this  particular  problem.  Therefore,  a  data  base 
was  generated  with  the  program  and  the  results  modeled  for  use  in  the 
TORCH  program.  The  most  sensitive  thermodynamic  property  for  this 
problem  is  the  flame  temperature.  It  is  assumed  that  this  temperature 
can  be  expressed  explicitly  as  a  function  of  quality  and  velocity  and 
implicitly  as  a  function  of  mass  ratio.  The  velocity  as  used  here  is 
a  measure  of  the  total  enthalpy  of  the  propane  jet.  The  flame  tempera¬ 
ture  may  then  be  expressed  as: 


T  =  T  (X,  U;  m) 


A  Taylor  series  expansion 


T  (X,  U;  m)  =  T  (0,  0;  m)  + 


+ 


(8) 


stipulates  that  the  flame  temperature  for  any  quality  or  velocity  may  be 
found  from  the  temperature  and  its  derivatives  evaluated  at  a  quality  of 
zero,  pure  liquid,  and  zero  velocity.  The  reference  temperature  and  its 
derivatives  in  this  case  remain  a  function  of  mass  ratio.  The  combustion 
program  was  exercised  with  a  variety  of  mass  ratios,  qualities,  and  velo¬ 
cities.  It  was  found  that  for  any  given  mass  ratio,  the  change  in  tem¬ 
perature  for  any  given  change  in  quality  or  velocity  was  linear,  indica¬ 
ting  that  the  first  order  Taylor  series  to  be  a  good  approximation. 
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Likewise,  the  enthalpy  can  be  written  as: 


h  (X,m  )  =  h  (0;  m)  + 


X 


(9) 


Molecular  weight,  MW,  and  specific  heat  ratio,  y,  were  found  to  depend 
on  the  mass  ratio  but  to  be  virtually  independent  of  quality.  The 
remaining  properties  such  as  density  and  sound  speed  can  be  computed 
assuming  the  ideal  gas  law  since  even  pure  propane  vapor  at  atmospheric 
pressure  and  temperature  will  behave  in  this  manner. 


The  functional  dependence  of  temperature,  enthalpy,  molecular 
weight,  and  specific  heat  on  mass  ratio  were  computed  for  mass  ratios 
less  than  stoichiometric  (m  =  17.3)  and  also  for  fuel  lean  mixtures.  It 
was  found  that  simple  polynomial  expressions  could  be  used  for  the  entire 
range  of  mass  ratios  with  the  exception  of  the  ratio  of  specific  heats, 
which  had  to  be  fitted  in  two  regions.  The  results  are  given  in  Table  1 


Table  1.  Polynomial  Coefficients,  a^ 
4 


n  =  0 


n 

T  (0,  0 

;  m) 

MW  (m) 
(g/mole) 

Y  (m) 

(m  <  17.3) 

Y  (m) 

(m  >  17.3) 

h  (0;  m) 
(cal/ g) 

0 

-9.3681 

E+2 

1.283  E+l 

1.3635 

1.3296  E-l 

-1.939  E+2 

1 

4.3158 

E+2 

1.8229  E+0 

-8.254  E-3 

1.2014  E-l 

2.2579  E+l 

2 

-1.9206 

E+l 

-6.9251  E-2 

0 

-4.2735  E-3 

-1.2611  E+0 

3 

2.5115 

E-l 

8.6095  E-4 

0 

5.1369  E-5 

3.3207  E-2 

4 

0 

0 

0 

0 

5.303  E-4 

The  functional  form  of  the  partial  derivations  in  Equations  (8)  and  (9) 
were  found  to  be: 


-2.4734  m  +  50.73  (°K) 


(10) 
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(11) 


(-1.41  m  +  28.894)  x  10'5  (°k  sec/cm) 


=  0.182  h  (0;  m)  (cal/g) 


(12) 


The  coefficients  of  the  various  terms  were  determined  by  minimizing  the 
square  of  the  deviation  of  the  assumed  equation  from  the  numerical  re¬ 
sults  of  the  combustion  program.  A  comparison  between  the  thermodynamic 
quantities  computed  by  the  NASA  program  and  the  fitted  polynomials  is 
shown  in  Figures  4-7. 


D.  Turbulent  Mixing 

The  mass  ratio,  m,  is  obviously  an  important  parameter.  The  amount 
of  air  entrained  is  a  function  of  downstream  position  and  jet  velocity. 
Analysis  of  the  centerline  velocity  decay  and  lateral  spreading  of  the 
jet  follows  the  basic  theory  of  Schlichting^as  modified  by  more  recent 
work. 


The  centerline  velocity,  Uc,  of  the  jet  is  reduced  by  the  entrain¬ 
ment  of  air  due  to  turbulent  mixing.  The  decay  may  be  expressed  quanti¬ 
tatively  as: 


U 


c 


1  -  exp 


(13) 


Here,  K  is  a  constant  related  to  the  eddy  viscosity;  X  is  the  nondimen- 
sional  axial  coordinate  expressed  in  jet  radii;  p  is  the  ratio  of  the 
ambient  air  density  to  the  jet  exit  density;  and  1.  is  a  jiondimensional 

potential  core  length  which  has  the  universal  value  o^  0  79  Tie  nara-i 
eter  K  varies  with  jet  Mach  number,  M  ,  and  is  assumed  to'  be*  coristahrt. 
The  values  used  in  this  analysis  are  according  to  Witze*  and  are  given 
by: 


H.  Sahlichting,  Boundary  Layer  Theory,  6th  Edition,  Chapter  24,  McGraw- 
Hill  Book  Company,  1968. 

12 

Colemans  DuP.  Donaldson  and  K.  Evan  Grey ,  "Theoretical  and  Experimental 
Investigation  of  the  Compressible  Free  Mixing  of  Tijo  Dissimilar  Gases," 
AT A A  Journal,  Vol.  4,  No.  11,  November  1966,  pp.  2017-2026 . 
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K  =  0.063 


-0.15 

1  M  >  1 

o 


and 


(14) 


K  =  0.08  (1  -  0.16  M  )  p  ~0-22  M  <  1  (15) 

o  o  — 

The  jet  is  composed  of  a  core  region  where  the  velocity  is  constant. 
Beyond  the  length  of  this  core,  the  centerline  velocity  is  given  by  Equa¬ 
tion  (13).  Away  from  the  centerline,  i.e.,  in  the  radial  direction,  the 
velocity  in  the  axial  direction  decays  exponentially: 


U  (r)  =  Uc  exp 


(16) 


where  r  ^  is  the  radial  position  where  the  velocity  is  one-half  of  the 
centerline  velocity. 

At  any  downstream  location,  the  jet  w.yi  spread  as  the  centerline 
velocity  decays.  According  t0  Kleinstein  the  radius  of  the  jet  to  the 
point  where  the  velocity  is  one-half  of  Uc  is: 


-  -1/2 

(X)  =  0.074  X  -  0.70  (p) 


(17) 


In  addition,  the  total  enthalpy  will  decay  a  form  similar  to  the  cen¬ 
terline  velocity  and  is  given  by  Kleinstein  as: 


Hc  (X) 

fi  - 

o 


exp 


1 


(0.102  /p  X  -0.7) 


(18) 


From  liquation  (18),  the  temperature  decay  along  the  centerline  may  be 
computed  from  the  energy  equation 


C  T  fX)  =  II  ( X) 
pc  c 


(19) 
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All  of  the  estimates  made  in  this  section  are  based  on  the  free  jet  expan¬ 
sion  by  turbulent  mixing.  The  presence  of  a  target  downstream  may  influ¬ 
ence  measurements  made  in  the  vicinity  of  the  target.  Therefore,  it  is 
necessary  to  make  an  estimate  of  the  wall  influence. 

It  was  shown  by  Rand"  that  the  centerline  velocity  in  the  presence 
of  a  wall  may  be  estimated  by 


U 

Cwall 


«  ( 

cno  wall\ 


1  -  exp 


(20) 


where  Z  is  the  distance  from  the  wall  and  rQ  is  the  jet  radius  (r  5  is 
a  convenient  value  to  use)  if  no  wall  had  been  present.  Sample  calcula¬ 
tions  indicate  that  measurements  by  probes  30  cm  (one  foot)  from  the  wall 
will  be  influenced  less  than  5  perceni  by  the  wall  for  fuel  rich  mixtures. 
However,  for  very  lean  mixtures,  the  radius  of  the  combustion  products  is 
sufficiently  large  for  the  ratio  (Z/r0)  used  in  Equation  (20)  to  become 
quite  small.  Under  these  circumstances,  the  influence  of  the  wall  is 
more  significant  and  the  core  length  is  sufficient  to  engulf  the  probe. 

The  exponential  form  of  Equation  (16)  is  such  that  the  maxium  velo¬ 
city  occurs  at  only  one  point  and  decays  to  one  half  the  maximum  value  at 
r  5  It  would  appear  reasonable  to  expect  the  measured  velocities  would 
be  more  related  to  the  average  value  than  the  local  maximum.  Therefore, 
average  velocity,  0,  is  defined  over  that  interval  where  the  velocity  is 
equal  to  or  greater  than  one  half  of  the  centerline  velocity. 


U  = 
c 


.5 


r  U  (r)  dr 


(21) 


Utilizing  Equation  (16),  the  average  velocity  is  given  by  the  integral 

1 

* 

0  =  Uc  2J y  exp  [ - ( «-n 2)  y2]  dy  (22) 

o 


so  that 


0  =  .72  0£  (23) 

This  simple  result  will  be  used  later  during  comparisons  with  experimentally 
obtained  velocities  and  heat  fluxes. 


31 


E.  Heat  Transfer 


The  purpose  of  the  facility  is  to  test  various  insulation  systems 
suitable  for  use  on  railroad  cars.  In  order  to  complete  the  analysis  of 
this  test  configuration,  the  heat  transfer  to  the  target  must  be  computed. 
Measurement  of  heat  transfer  is  accomplished  with  a  commercially  avail¬ 
able  probe!  however,  some  of  the  coefficients  are  dependent  on  geometry 
and  flow  conditions.  Therefore,  a  brief  discussion  of  the  heat  transfer 
to  the  target  plate  is  warranted. 

In  choosing  an  appropriate  analytic  expression  for  the  heat-transfer 
coefficient,  a  number  of  physical  processes  must  be  considered.  Imme¬ 
diately  upon  leaving  the  orifice,  and  after  expansion  to  atmospheric 
pressure,  the  jet  begins  to  entrain  air.  This  entrainment  leads  to  tur¬ 
bulence,  which  is  further  enhanced  by  the  combustion  process.  Increasing 
free  stream  turbulence  has  a  direct  effect  on  local  heat  transfer  coeffi¬ 
cients.  However,  for  both  fully  developed  turbulent  jets  and  for  initially 
laminar  jets  that  become  turbulent  as  a  result  of  jet  mixing  (e.g., 
entrainment  of  air) ,  the  intensity  of  turbulence  appears  to  be  uniquely 
determined^  by  the  jet  Reynolds  number  and  the  dimensionless  jet  length, 
expressed  in  terms  of  effective  orifice  diameters.  Thus,  the  data  on 
heat-transfer  rates  can  be  correlated  very  effectively  without  any  sep¬ 
arate  parameter  to  characterize  turbulence. 

The  jet  hits  the  plate  and  a  stagnation  boundary  layer  is  formed. 
Figure  8.  In  the  stagnation  region,  the  flow  field  is  laminar  and  the 
boundary  layer  is  thin,  j  •>  Since  the  boundary  layer  is  thin,  the 
heat-transfer  coefficient  can  be  large.  It  has  been  noticed  that  "... 
the  disappearance  of  the  pressure  gradients  which  exist  in  the  vicinity 
of  the  stagnation  point  ...  serve(s)  to  stabilize  the  laminar  boundary 
layer,  in  spite  of  locally  already  high  turbulence  levels  in  the  free 
stream.  Thus,  it  is  only  at  the  outer  edge  of  the  jet  deflection  region 
that  conditions  are  conducive  for  a  ti^psition  from  a  laminar  to  a  tur¬ 
bulent  boundary  layer  to  take  place." 
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Figure  &.  The  Impinging  Jet 


In  the  present  analytic  model  of  atcombusting  propane  jet  interact¬ 
ing  with  a  flat  plate,  the  reaction  zone  is  assumed  to  be  physically  sepa¬ 
rate  from  the  jet/plate  interaction.  In  reality,  this  may  not  be  true  and, 
therefore,  an  expression  for  chemical  reactions  might  need  to  be  included 
in  the  computation  of  tl^  l^^at-transfer  coefficient.  However,  the  data 
from  some  investigators  have  shown  that  good  agreement  can  be  obtain¬ 

ed  with  experimental  data  even  if  the  chemical  kinetics  may  make  the  ana¬ 
lytic  prediction  a  worse  predictor  of  the  heat-transfer  coefficient  (except 
perhaps  for  a  catalytic  surface).  These  findings  may  be  explained  by  the 
fact  that  the  laminar  boundary  layer  is  thin;  thus,  the  number  of  reactions 
contributing  to  the  heat  transfer  in  the  boundary  layer  is  very  small  and, 
therefore,  can  be  neglected  when  computing  the  heat-transfer  coefficient. 

Thus,  an  analytic  expression  for  the  Nusselt  number,  Nu,  in  terms  of 
the  Prandtl  number,  Pr,  and  Reynolds  number,  Re,  can  be  derived  from  lam¬ 
inar,  stagnation  flow: 

Nu  =  Constant  •  Pr^'4  Re  (24) 

K  K 

where 


Nu  =  hR/k  ( 2S) 


Here,  h  is  the  convective  heat  transfer  coefficient,  R  is  the  stagnation 
radius,  and  k  is  the  thermal  conductivity  of  the  boundary  layer.  The 
Reynolds  number  is  computed  from  free  stream  variables,  while  the  Prandtl 
number  and  thermal  conductivity  are  computed  using  the  average  temperature 
of  the  free  stream  and  wall  temperatures.  The  constant  of  proportionality 
is  dependent  upon  the  geometry;  for  a  three-dimensional  jet  of  radius  R 
impi^g;ng  upon  a  flat  plate,  the  constant  has  a  value  of  the  order  of 
0.5  J  .  The  stagnation  radius,  R,  is  taken  to  be  the  width  of  the  jet  at 
half-velocity,  Equation  (IT),  and  the  velocity  used  in  computing  the 
Reynolds  number  is  the  average  value  defined  by  Equation  (23).  The  heat- 
transfer  coefficient  can  then  be  written: 


h  =  0.45  k  (Uc/RvS'S  Pr°-4  (26) 

where  v  denotes  the  kinematic  viscosity  of  the  fluid. 

Since  the  stagnation  radius  and  average  velocity  are  dependent  upon 
the  orifice-to  plate  distance,  it  can  be  seen  that  the  heat -transfer  to 
the  plate  is  also  dependent  on  this  distance.  As  the  centerline  velocity 
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decreases  exponentially  with  distance,  and  R  increases  linearly  with  dis¬ 
tance,  the  heat  transfer  to  the  plate  is  obviously  reduced  as  the  plate 
is  moved  away  from  the  orifice.  From  Equation  (26),  it  is  obvious  that 
the  properties  of  the  combusting  jet  such  as  gas  velocity,  gas  density, 
viscosity,  flame  conductivity  and  flame  temperature  (the  fluid  properties 
are  temperature  dependent  and  the  heat  transfer  depends  upon  the  differ¬ 
ence  in  flame  temperature  and  plate  temperature)  are  important  in  comput¬ 
ing  the  actual  heat  transfer  from  the  jet  to  the  plate. 

Modeling  assumptions,  uncertainties  in  calculating  some  of  the  jet 
parameters,  and  difficulties  in  taking  experimental  measurements  of  the 
"torch"  require  that  Equation  (26)  be  altered  to  determine  the  constant 
of  proportionality  which  allows  the  best  correlation  with  the  experimental 
data. 


Outside  the  stagnation  region,  the  convective  heat  transfer  to  the 
plate  changes  form  due  to  the  transition  from  laminar  to  turbulent  flow. 
Thus,  for  the  wall  jet  region  depicted  in  Figure  8,  the  appropriate  analy¬ 
tic  expression  for  the  convective  heat  transfer  is  for  a  turbulent  bound¬ 
ary  layer  over  a  flat  plate:  ° 


Nur  =  0.0288  ReR4/5  Pr1/3,  r  >  R 


(27) 


where,  now,  the  Nusselt  and  Reynolds  numbers  are  a  function  of  the  radial 
distance  from  the  stagnation  point.  The  wall  jet  velocity  decays  accord¬ 
ing  to: 


Ur  =  U  ^  1  -  exp 


wir 


(28) 


Hence,  the  heat  transfer  coefficient  can  be  computed: 


h  =  0.0243 


Pr 


1/3 


1  -  exp  I- 


mr 


r  >  R  (29) 


Again,  due  to  uncertainties  in  various  parameters,  the  constant  of  propor¬ 
tionality  is  altered  for  the  best  correlation  with  experimental  data.  The 
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heat  flux  to  the  plate,  q,  consists  of  a  convective  and  radiative  compo¬ 
nent,  and  can  be  written  as: 


q  -  hrCTp  -  Tp) 


+  oe  e 
P  F 


crt 


tp4) 


(30) 


Here,  the  subscripts  F  and  P  refer  respectively  to  the  flame  and  plate; 
e  are  their  respective  emissivities;  a  is  the  Stefan-Boltzmann  constant. 
The  radiation  exchange  between  the  flame  and  the  plate  is  approximated 
by  the  radiation  exchange  between  two  infinite  plates.  The  flame  emis- 
sivity  is  a  function  of  quality  -  for  a  predominately  liquid  flame  (qual¬ 
ity  'v  0),  Cp  approaches  1.0  because  the  fuel  rich  propane/air  mixture 
produces  considerable  radiation.  9  For  the  diffusion  flames,  it  has  been 
noted  that  the  tendency  for  carbon  and  soot  formation  increases  with  mass 
flow  and  for  very  fuel  rich  flames  which  results  in  a  larger  radiative 
component.  The  flame  is  optically  thinner  for  a  "vaporous"  flame  and 
the  radiative  component  to  the  heat  transfer  is  substantially  reduced. 

In  summary,  the  heat  transfer  from  the  burning  propane  torch  to  the 
plate  consists  of  a  radiative  and  convective  component,  Equation  (30), 
where  the  convective  component  is  given  by: 


h 

r 


k  0^04 

h  =  a  4  (U  R/v)  Pr  r  <  R 

r  R  c  —  . 


n  k  ...  D/  ,..0.8  0.33  .  . 

8  (U^R/v)  Pr  |  1  -  exp 


v0-4 

[2  (*-)]! 


r  >  R 


(31) 

(32) 


where  a  and  8  are  found  from  correlation  with  the  experimental  data. 

In  the  experiments  performed  by  the  BRL,  a  calorimeter  manufactured 
by  Thermogage  was  placed  approximately  along  the  centerline  of  the  torch. 
The  sensing  element  is  in  the  tip  of  a  cone  where  the  apex  of  the  cone  has 
been  "rounded,"  and  has  a  radius  of  approximately  0.32  cm.  In  order  to 
compare  the  experimental  data  with  the  theory,  the  heat-transfer  coeffi¬ 
cient  of  the  torch  to  the  Thermogage  probe  must  be  estimated.  Equation 
(31)  is  applicable  but  the  constant  must  now  be  evaluated  for  the  geometry 
of  the  probe.  The  result  of  this  analysis  yields: 


A.  G.  Gay don  and  H.  G.  Wolfhard ,  Flames,  Their  Structure ,  Radiation  and 
Temperature ,  4th  Edition,  "A  Halsted  Press  Book,  "  John  Wiley  <$  Sons,  1979. 


h  .  .  =  1.07  (U  /R  .  V)0'5  P  0,4 

torch+probe  c  probe  r 


(33) 


where  R  .  =  0.32  cm. 

probe 


F.  Parametric  Analysis 

The  analysis  procedures  described  here  have  been  programmed  for  the 
purpose  of  rapid  numerical  evaluation.^  The  results  shown  in  Figure  2  are 
typical  of  the  flow  properties  at  the  exit  of  the  orifice.  It  may  be  seen 
that  the  jet  exit  velocity  is  a  strong  function  of  reservoir  quality  and 
pressure. 

The  results  of  the  combustion  analysis  indicated  that  the  most  sig¬ 
nificant  thermodynamic  property  is  flame  temperature.  This  property  is 
most  sensitive  to  the  mass  ratio  of  air  to  fuel  and  relatively  insensitive 
to  other  parameters  such  as  reservoir  pressure  or  quality.*  Since  the 
target  is  engulfed  in  the  jet  formed  by  the  combustion  products,  and 
since  these  properties  are  independent  of  quality,  the  computed  flow  con¬ 
ditions  (per  our  modeling  assumptions)  at  the  target  are  also  independent 
of  quality. 

A  parametric  study  was  conducted  for  the  mass  ratios  and  reservoir 
pressures  of  interest  assuming  a  quality  of  one  (pure  vapor).  Although 
many  flow  properties  are  computed,  the  temperature  and  the  heat  flux  at 
the  probe  in  front  of  the  target  are  typical.  These  results  are  pre¬ 
sented  in  Figures  9  and  10  and  graphically  demonstates  the  insensitivity 
of  these  results  to  reservoir  pressure ,**The  points  on  these  figures  rep¬ 
resent  pressures  between  0.827  and  1.65  MPa  (120  and  240  psia)  which  spans 
the  region  of  interest  at  the  torch  simulation  facility. 

Another  result  of  this  parametric  study  was  the  observation  that  for 
large  mass  ratios,  the  core  length  of  the  jet  formed  by  the  combustion  pro¬ 
ducts  was  sufficiently  long  to  engulf  the  measuring  probe.  This  effect  is 
observed  for  fuel  lean  mixtures  and  is  caused  by  the  addition  of  large  quan 
tities  of  air  which  slows  the  original  propane  jet  and  increases  the  radius 
of  the  mixture  to  the  point  that  the  core  length  predicted  by  diffusion 
theory  exceeds  the  distance  to  the  target.  Therefore,  the  reduced  tempera¬ 
tures  shown  in  Figure  10  are  more  indicative  of  an  excess  of  air  during  com 
bustion  than  cooling  during  the  diffusion  process. 

— - 

This  model  assumes  a  simplified  mixing  theory  where  the  fuel  and  air  are 
thoroughly  mixed  before  combustion.  With  this  assumption ,  the  flame  tem¬ 
perature  was  found  to  be  essentially  independent  of  quality.  A  more  de¬ 
tailed,  two-dimensional  analysis  would  be  required  to  account  for  the  two- 
phase  mixture  across  the  width  of  the  jet. 

Note  that  for  air-to-fuel  ratios  greater  than  stoichiometric,  that  the 
temperature  is  virtually  independent  of  pressure,  and  is  relatively  in¬ 
sensitive  to  pressure  for  ratios  lower  than  stoichiometric . 
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Figure  9.  Effect  of  Mass  Ratio  on  Heat  Flux  for  Various  Pressures 
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Ill.  DATA  ANALYSIS 


A  series  of  tests  were  conducted  at  the  torch  simulation  facility  in 
Pueblo,  Colorado,  from  28  October  to  2  November  1978.  The  configuration 
was  as  depicted  in  Figure  1.  Valves  A  and  B  were  electrically  driven,  man¬ 
ually  controlled  gate  valves  that  were  intended  to  control  the  quality  of 
the  effective  reservoir  upstream  of  the  orifice.  The  target  was  instru¬ 
mented  with  a  number  of  Tempellets  positioned  just  in  front  of  a  steel 
plate;  a  Hastings  velocimeter  probe  and  a  Thermogage  calorimeter  were 
centered  30  cm  in  front  of  the  plate.  Analog  data  were  taken  as  a  function 
of  time,  converted  to  digital  signals  and  recorded.  Data  were  taken  in¬ 
dicating  both  valve  rotation  positions,  pressure  upstream  of  the  orifice, 
velocity  and  heat  flux  at  the  probe  position.  The  recorded  data  from  13 
tests  have  been  collected  and  are  contained  in  Appendix  A  of  this  report. 

The  purpose  of  this  study  is  to  compare  the  analytic  model  predic¬ 
tions  to  measured  experimental  quantities.  The  data  contained  in  Appendix 
A  is  assumed  to  be  accurately  measured  and  calibrated.  An  attempt  has  been 
made  to  estimate  the  uncertainty  of  the  heat  flux  measurements  due  to  the 
fluctuating  nature  of  the  torch.  This  estimate  was  made  by  analyzing  the 
high  frequency  variability  of  the  data. 

A.  Valve  Data 

Two  identical  valves  were  used  to  control  the  flow  of  propane  to  the 
orifice.  Valve  A  controlled  the  liquid  phase  and  valve  B  the  vapor  phase. 
These  valves  have  a  two  inch  diameter  orifice  and  are  positioned  in  a  four 
inch  line.  It  was  recognized  that  some  rotation  of  the  valve  was  necessary 
before  any  flow  could  occur  due  to  the  design  of  the  valve.  It  was  also 
recognized  that  as  the  valve  was  rotated  and  the  flow  increased,  a  point 
would  be  reached  where  the  small  nozzle  at  the  exit  of  the  manifold  would 
choke  the  flow  and  allow  the  measured  pressure  to  approach  the  tank  pres¬ 
sure.  Further  rotation  of  the  valve  would  not  influence  the  measured  pres¬ 
sure  or  any  other  flow  parameter. 

No  calibration  data  was  supplied  with  the  experimental  records  so 
that  there  is  no  way  to  directly  compute  the  quality  of  the  reservoir  from 
the  recorded  valve  positions.  However,  a  study  of  the  recoids  in  Appendix 
A  revealed  that  a  number  of  tests  were  conducted  with  the  liquid  valve,  A, 
in  the  completely  closed  position.  These  four  tests,  147,  148,  149,  and 
150,  were  examined  in  detail  to  obtain  the  valve  opening  characteristics. 
All  data  is  scaled  to  values  between  zero  and  one  thousand  for  convenience 
in  presentation.  In  the  case  of  a  valve,  this  scale  corresponds  to  rota¬ 
tion  angles  from  zero  to  90°.  The  pressure  is  measured  just  upstream  of 
the  exit  nozzle  which  has  an  orifice  diameter  of  0.95  cm.  The  controlling 
valves  have  a  maximum  orifice  size  of  approximately  5  cm  when  fully  opened. 
When  the  calve  is  sufficiently  open,  the  measured  pressure  will  be  the  tank 
supply  pressure  and  the  exit  nozzle  will  serve  to  limit  the  flow. 

In  order  to  assess  the  effectiveness  of  valve  rotation,  the  ratio  of 
the  measured  pressure  upstream  of  the  exit  nozzle  and  the  tank  supply 
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pressure  is  shown  in  Figure  11  as  a  function  of  measured  valve  rotation. 

The  data  was  taken  from  tests  147  through  150  during  which  the  liquid  valve 
was  completely  closed.  The  pressure  response  of  Valve  B  is  surprisingly 
linear  with  valve  rotation.  Of  particular  interest,  at  positions  below 
176,  the  valve  may  be  considered  completely  closed.  Also,  at  positions 
above  410,  the  valve  is  open  sufficiently  for  the  tank  pressure  to  be 
measured  at  the  exit  nozzle.  This  observation  is  significant  in  that  cer¬ 
tain  conclusions  can  be  drawn  about  the  quality  of  the  reservoir  conditions. 

Valve  A  which  controls  the  flow  of  liquid  propane  to  the  nozzle  exit 
is  identical  to  valve  B.  Due  to  the  difference  in  densities  between  liquid 
and  vapor  at  high  pressures,  the  onset  of  flow  through  valve  A  should  re¬ 
quire  the  same  or  slightly  larger  rotation  than  valve  B.  Therefore,  at 
valve  settings  below  176,  no  liquid  is  being  supplied  to  the  manifold  and 
the  resulting  reservoir  upstream  of  the  nozzle  exit  is  supplied  by  the  sat¬ 
urated  vapor  supply  only. 

An  examination  of  the  experimental  records  contained  in  Appendix  A 
indicates  that  in  the  majority  of  cases,  the  liquid  valve  may  be  considered 
to  be  completely  closed  and  the  vapor  valve  completely  open.  Only  tests  141 
and  146  indicate  that  valve  A  was  open  sufficiently  to  provide  any  flow 
of  liquid.  Therefore,  the  analysis  used  for  comparison  with  experimental 
data  have  an  assumed  quality  of  one. 

B.  Velocity  and  Calorimeter  Data 


At  this  point,  it  is  appropriate  to  review  a  typical  experimental 
record  to  understand  the  manner  in  which  the  data  will  be  examined.  Test 
137  is  typical  and  will  be  discussed  in  detail.  The  relationship  between 
the  various  symbols  shown  and  the  values  plotted  is  given  at  the  beginning 
of  the  Appendix.  The  five  quantities  measured  were:  A,  the  liquid  valve 
opening;  B,  the  vapor  valve  opening;  G,  orifice  pressure  (psig) ;  N,  the 
heat  flux  (Btu/ft^  sec);  and  0,  the  velocity  (ft/min  at  standard  tempera¬ 
ture  and  pressure) . 


The  initial  conditions  indicated  on  this  record  are  as  follows: 


Quantity  Value 


Remarks 


A 

B 

G 


N 


0 


125  Rotated  11.25°  -  no  flow 
410  Rotated  36.9°  -  fully  open 

0  Indicates  no  fuel  at  ori¬ 

fice;  valve  C  (Figure  1) 
probably  closed 

0  No  heat  transfer  indicates 

either  no  flow  or  no  igni¬ 
tion  of  mixture 

150  Velocity  measured  is  300 

ft/min  (2.96  knots).  This 
corresponds  to  the  observed 
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wind  conditions  at  the  time 
in  the  direction  of  the  tar¬ 
get 

At  1.2  minutes  into  the  test,  valve  C  is  opened  as  evidenced  by  the 
sudden  increase  in  orifice  pressure.  Significant  increase  in  measured 
velocity  is  indicated  but  without  heat  transfer.  Flow  is  obviously  an 
unignited  mixture  of  propane  and  air. 

At  3.5  minutes  into  the  test,  valve  C  is  closed  and  at  4  minutes  the 
data  recorder  is  interrupted  and  then  restarted  at  9.5  minutes.  At  10.1 
minutes,  valve  C  is  again  opened  and  mixture  ignited.  Both  velocity  and 
heat  flux  are  shown  to  be  fluctuating  with  peak  values  of  1540  ft/min 
and  43.5  Btu/ft^  -sec  respectively.  The  test  is  terminated  after  14.5 
minutes  by  closing  valve  C. 

Each  record  must  be  examined  individually  to  extract  meaningful  in¬ 
formation.  The  velocity  measuring  device  has  a  relatively  slow  response 
time  so  that  fluctuations  are  generally  considered  to  be  an  integrated  ef¬ 
fect  of  locally  turbulent  conditions  caused  by  a  variety  of  sources  in¬ 
cluding  gusting  winds.  The  peak  values  are  of  most  significance  and  are 
listed  in  Table  2.  In  some  tests,  the  flow  conditions  were  changed  during 
the  course  of  the  test.  Therefore,  multiple  data  points  are  listed  in 
those  cases.  The  calorimeter  data  shown  in  the  experimental  records  ap¬ 
pears  to  be  the  RMS  value  of  the  analog  data..20  The  fluctuations  in  the 
analog  data  are  of  a  relatively  high  frequency  with  a  peak  to  peak  ampli¬ 
tude  on  the  order  of  20  percent  of  the  mean  value.  Both  velocity  and  heat 
flux  data  have  been  converted  to  cgs  units  so  that  direct  comparison  with 
the  analysis  contained  in  Section  II  is  possible.  These  dimensional  values 
are  also  listed  in  Table  2. 

IV.  RESULTS  AND  DISCUSSION 

In  order  to  predict  the  velocity  and  heat  flux  at  the  target  using 
the  analysis  procedures  contained  in  Section  II,  it  is  necessary  to  have 
some  knowledge  of  the  mass  ratio  of  air  to  propane.  This  parameter  domi¬ 
nates  the  combustion  process  which  in  turn  controls  the  downstream  flow 
properties.  The  theoretical  results  already  presented  as  Figures  9  and  10 
are  replotted  as  a  function  of  reservoir  pressure  for  various  mass  ratios 
in  Figures  12  and  13.  In  this  form  the  experimental  values  contained  in 
Table  2  may  be  presented  since  both  ordinate  and  abscissa  are  measured 
quantities . 

The  velocity  data  contained  in  Table  2  is  superimposed  on  the  anal¬ 
ysis  results  in  Figure  12.  Different  symbols  have  been  used  to  identify 
those  data  obtained  with  the  liquid  valve  completely  closed.  It  can  be 
seen  from  this  figure  that  all  usable  data  falls  within  a  relatively 

20Charles  Anderson ,  private  communications. 
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Figure  12.  Effect  of  Reservoir  Pressure  on  Average  Velocity  at  Plate 
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Figure  13.  Effect  of  Reservoir  Pressure  on  Heat  Flux  to  Probe 


narrow  band  bounded  by  mass  ratios  from  17.5  to  24.5.  At  the  inception 
of  this  program,  it  was  expected  that  combustion  would  occur  first  with 
a  fuel  rich  mixture  which  is  represented  by  mass  ratios  less  than  17.3. 
However,  the  data  suggests  that  combustion  occurs  in  this  facility  only 
with  a  lean  mixture.  In  practice,  the  combustion  of  a  fuel  jet  is  actually 
a  two-dimensional  phenomenon.  Ignition  probably  occurs  at  the  periphery 
of  the  jet  where  lean  mixture  exists.  Once  ignited,  the  flame  may  pro¬ 
pagate  to  the  centerline  so  long  as  the  mixture  is  sufficient  to  support 
combustion.  The  one-dimensional  measurement  of  flow  properties  is  actually 
a  measurement  of  the  integral  effect  over  the  width  of  the  jet.  For  lean 
mixtures,  the  assumption  of  thorough  air-fuel  mixing  before  combustion 
appears  to  be  a  relatively  valid  approximation  (i.e.,  comparison  of  analytic 
predictions  with  the  experimental  data).  Thus,  the  validity  of  the  com¬ 
plete  mixing-before-combustion  assumption,  along  with  the  integral  effect, 
allow  a  one-dimensional  fuel  lean  mixture  model  to  yield  realistic  predic¬ 
tions  . 

Since  the  data  obtained  from  those  tests  with  the  liquid  valve  com¬ 
pletely  closed  does  not  overlap  the  other  data,  it  may  be  concluded  that 
a  high  quality  flow  such  as t  hat  produced  by  a  saturated  vapor  source  will 
bum  closer  to  stoichiometric  condition  than  more  fluid  flames.  The  bal¬ 
ance  of  the  data  suggests  that  the  addition  of  very  little  liquid  will  in¬ 
crease  the  mass  of  air  required  for  combustion  disproportionately.  The 
theoretical  observation  that  the  resulting  flow  parameters  are  relatively 
independent  of  reservoir  pressure  is  confirmed  by  the  experimental  obser¬ 
vations. 

The  variability  in  the  peak  velocity  data  is  well  within  the  expected 
range  of  scatter  due  to  wind  gusts  and  the  difficulty  in  aiming  the  propane 
source.  Bouyancy  effect  of  the  flame  coupled  with  less  than  ideal  wind  con¬ 
ditions  make  it  very  difficult  to  center  the  point  of  impingmeht  within  a 
circle  of  radius  8  cm  at  a  distance  of  365  cm.  This  is  the  radius  of  the 
theoretical  half  velocity  point  under  stoichiometric  conditions.  Unfortun¬ 
ately,  the  variability  of  this  data  is  sufficient  to  mask  any  subtle  effects 
of  pressure.  At  the  same  time,  the  data  indicates  that  the  mass  ratio  is 
influenced  by  very  small  changes  in  the  liquid  valve  setting. 

The  calorimeter  data  contained  in  Table  2  are  superimposed  on  the  ana¬ 
lytic  results  in  Figure  13.  In  addition,  an  error  bar  has  been  included 
which  is  ±  10  percent  of  the  measured  flux.  This  is  based  on  the  fact  that 
analog  data  indicates  an  oscillation  of  this  order  of  magnitude.  In  this 
figure,  only  theoretical  results  for  fuel  lean  mixtures  are  presented. 

These  results  are  dependent  of  the  radius  of  the  probe  which  was  estimated 
to  be  0.32  cm.  Theoretical  agreement  could  be  obtained  by  increasing  this 
number  slightly;  however,  in  the  absence  of  details  of  the  probe,  these  re¬ 
sults  are  considered  adequate  to  explain  the  heat  transfer  mechanism.  The 
experimental  data  does  not  seem  to  support  the  reduction  in  heat  flux  for 
increased  mass  ratios.  This  is  believed  to  be  due  to  the  increase  in 
emissivity  with  increasing  amounts  of  fluid.  This  feature  has  been  repor¬ 
ted  by  observers  but  has  not  been  taken  into  account  in  the  analysis. 
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As  in  the  case  of  the  velocity  measurements,  a  certain  amount  of  varia¬ 
bility  in  the  data  may  be  attributed  to  the  aim  of  the  propane  jet  and 
local  wind  conditions.  In  any  event,  the  experimental  data  suggest  that 
the  analytic  model  can  predict  the  incident  velocity,  heat  flux  and  tem¬ 
perature  of  a  reacting  propane  jet  with  a  quality  near  one. 

V.  RECOMMENDATIONS 

The  experimental  data  from  the  torch  simulation  facility  has  pro¬ 
vided  the  necessary  evidence  that  the  thermodynamic  and  fluid  dynamic 
models  developed  are  representative  of  the  tests.  The  analysis  may  be 
used  to  identify  where  the  experiment  may  be  improved  and  the  experimen¬ 
tal  data  suggests  where  certain  improvements  in  the  analysis  are  possible. 

The  torch  simulation  facility  may  be  improved  by  noting  the  sensitiv¬ 
ity  of  the  data  to  the  mass  ratio.  Although  the  combustion  process  is  rela 
tively  independent  of  pressure  and  quality,  the  mass  ratio  appears  to  be 
a  strong  function  of  the  liquid  injected  into  a  saturated  vapor  flow. 
Therefore,  better  metering  of  the  liquid  is  necessary  if  more  precise  in¬ 
formation  is  desired  from  this  test.  This  would  require  a  new  valve  or 
precision  orifices  to  be  installed  in  place  of  valve  A,  Figure  1. 

Much  of  the  variability  in  the  data  is  caused  by  the  inability  to  cen 
ter  the  jet  downstream  due  to  buoyance  and  wind  effects.  These  effects 
could  be  reduced  somewhat  by  conducting  bench  scale  testing  in  an  enclosed 
environment.  However,  there  is  no  evidence  to  indicate  that  scaling  of 
turbulent  diffusion  flames  is  feasible. 

The  model  has  assumed  thorough  mixing  of  the  propane  and  air  before 
combustion. Since  the  model  makes  predictions  which  are  in  agreement  with 
experimental  data,  this  assumption  appears  to  be  supported.  However,  the 
experimental  data  is  for  initial  qualities  near  1.0,  i.e.,  essentially 
vapor  flow  with  a  small  amount  of  liquid.  For  thorough  mixing,  the  two- 
phase  analysis  indicated  that  the  flame  temperature  is  virtually  in¬ 
dependent  of  quality.  Observations  of  the  test  clearly  indicate  that  for 
initial  qualities  near  zero  i.e.,  essentially  liquid  flow  that  the  com¬ 
bustion  process  is  definitely  two-dimensional;  it  is  a  fin  -lion  of  radial 
distance  from  the  jet  centerline  and  axial  distance  from  the  orifice.  A 
non-burning  core  of  propane  is  clearly  pvident  for  small  qualities. 

Additionally,  the  torch  radiates  significant  amounts  of  energy  when 
the  initial  quality  is  zero.  This  is  certainly  due  to  incomplete  com¬ 
bustion  processes  in  a  fuel  rich  zone  creating  carbon  particles  (soot) 
which  can  radiate  large  amounts  of  energy.  Indeed,  a  "liquid"  flame 
impacting  a  flat  plate  transfers  significantly  less  heat  to  the  plate 
compared  to  the  "vapor"  flame,  indicating  a  combination  of  a  cooler  ef¬ 
fective  flame  temperature  and  radiative  heat  losses  due  to  the  incomplete 
combustion  processes. 
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The  analysis  could  be  improved  by  allowing  the  emissivity  to  be  a 
function  of  mass  ratio.  This  could  be  accomplished  by  measuring  the  radia¬ 
tion  as  well  as  the  heat  flux  to  a  calorimeter.  In  the  present  analysis, 
the  emissivity  is  assumed  to  be  constant;  however,  as  more  fluid  is  added, 
the  increase  in  radiant  heat  is  quite  apparent  to  personnel  at  the  facility. 
The  heat  flux  due  to  convection  is  reduced  as  the  velocity  is  decreased 
with  increasing  fluid.  However,  if  the  radiant  energy  released  increases, 
the  net  measured  heat  flux  would  be  the  sum  of  the  two  components 
and  would  be  less  sensitive  to  increases  in  mass  ratio  for  small  changes 
of  quality  than  the  present  analysis  indicates. 

In  addition,  an  analysis  of  the  mixing  process  could  be  performed 
using  a  two-dimensional  code  (e.g.,  GENMIX^  j  to  examine  the  effects  of 
quality  on  the  mixing  process.  An  expected  result  would  be  that  the  effec¬ 
tive  flame  temperature,  and  consequently  the  heat  transfer,  would  be  found 
to  be  dependent  upon  the  quality. 

Further,  the  analysis  here  only  considered  one  location,  3.66m  (12.0 
ft),  downstream  from  the  orifice.  Data  at  other  locations  would  be  required 
to  ascertain  the  predictive  ability  of  the  present  model  at  other  downstream 
locations. 

In  conclusion,  it  has  been  shown  that  the  analysis  techniques  devel¬ 
oped  in  this  program  are  sufficiently  accurate  to  predict  the  observed  phe¬ 
nomenon  for  the  limited  range  of  an  initial  quality  near  one.  The  torch 
facility  is  generally  used  with  a  quality  near  one,  thus,  this  restriction 
is  not  too  severe.  To  fully  understand  the  torching  mechanism  for  various 
qualities,  additional  data  is  required  to  establish  the  mass  ratio  as  a 
function  of  liquid  valve  position;  and  an  analysis  of  the  aforementioned 
effects  of  mixing,  combustion,  and  emissivity  are  required.  For  the  data 
available,  the  mass  ratio  varies  from  17.3  (stoichiometric)  to  approximately 
21  when  the  liquid  valve  is  in  the  12  percent  position  and  the  vapor  valve 
is  in  the  40  percent  position.  There  is  no  reason  to  believe  that  this  re¬ 
lationship  is  linear  and  a  carefully  planned  series  of  tests  is  warranted 
to  provide  the  necessary  information  to  characterize  the  flow  condition. 
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APPENDIX  A 


RECORDED  DATA  FROM  TESTS 


NOMENCLATURE 


A:  Liquid  valve  opening 
B:  Vapor  valve  opening 
G:  Orifice  pressure  *  3  (PSIG) 

N:  Thermogage  calorimeter  *  10  (BTU/f t"2-sec) 

0:  Hastings  velocimeter  *  0.5  (ft/min:  air  @  STP) 
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APPENDIX  B 

TABULATION  OF  COMPUTER  MODEL 
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PROGRAM  HOTPLT 

DIMENSION  TS(100*9)*T( 100*8*  2 ) » HTC (100*8 ) 

DIMENSION  THEAD(125)*HTCBAR(10C*6) 

DIMENSION  P(20)*TK(20)*SF(20)*SG(20)*HF(  20)*  HG  ( 20 )*  VF(  20),  VG(  20) * 

♦  VEL  (  20  »  X(  20 ) » H ( 20  )*  RHO(  20) *  V  ( 20 )»  WDOTA(  20)*S(20) 

DIMENSION  PPR(10)»MQF(10)*VASTPC(10*1G)*QC(10*10) 

COMMON/ PRMTRS/TS»T» HTC *TRES*  TA 

REAL  MU9AR*KBAR»MUA*KI*KS*KSTL*MEFF»KAPPA*MW*MRATI0*MACHCP*MWCPC 

REAL  IGDST* MOF 

DATA  HTC*HTCBAR/1600*0«0/ 

DATA  T*TS, TINF/2401*549.1/ 

OPEN  (JNIT-l,MAME-'HaTIN.DAT»*TYPE«‘OLD* ) 

OPEN  (UNIT-2*NAME-'H0T0UT.0AT'*TYPE-*NEW» ) 

OPEN  (UNIT *4#  NAME* 'HOTLIP. DAT** TYPE- 'NEW  • ) 

C 

C  IN»UT  QUANTITIES  AND  DIMENSIONS 
C  I  THSTL-PLATE  THICKNESS  (FT) 

C  I  TA-TENPERATURE  OF  JET  TIP  (R) 

C  I  KI-THERMAL  CONDUCTIVITY  OF  INSULATION  ( BTU/ SE C-FT-F ) 

C  I  THI- 1 NSUL A  T I  ON  THICKNESS  (FT) 

C  I  KS-THERMAL  CONDUCTIVITY  OF  SOOT  ( BTU/ SEC-FT-F ) 

C  I  THS-SOOT  THICKNESS  (FT) 

C  I  TINF- AMBIENT  TEMPERATURE-INITIAL  PLATE  TEMPERATURE  (R) 

C  I  SETA-TILT  ANGLE  OF  PLATE  FROM  HORIZONTAL  (OEG) 

C  I  PLDST-D I  ST  ANC  £  FROM  NOZZLE  TO  PLATE  (FT) 

C  PBDS  T -D ISTANCE  FROM  NOZZLE  TO  PROBE  (FT) 

C  I  GDI  ST- DISTANCE  FROM  NOZZLE  TO  IGNITOR  (FT) 

C  I  Z-CQMPRESSI8ILITY  OF  PROPANE  AT  PEFF  AND  TA 

C  I  RT-NOZZLE  THROAT  RADIUS  (FT) 

C  I  PLATE-HALF  SPAN  OF  SMALLEST  PLATE  EDGE  (FT) 

C  I  DELT-TIME  INCREMENT  SELECTED,  MUST  Be  LESS  THAN  DELT1  AND 

C  DELT2  (SEC) 

C  I  E-EMISS IVITY  OF  SURFACE  MATERIAL 

C 

C  C  DELT1-L ARGEST  TI"E  INCREMENT  ALLOWED  TO  COMPUTE  T(I*J,L)  FOR 

C  THE  CENTER  ELEMENT  (SEC) 

C  C  DELT2 -L ARGEST  TIME  INCREMENT  ALLOWED  TO  COMPUTE  T(I,J*L)  FOR 

C  THE  GENERAL  ELEMENTS  (SEC) 

C  C  TRES-THERMAL  RESISTANCE  OF  SOOT  AND  INSULATION  (StC-FT/BTU) 

C  C  TS(I» J) -TEMPERATURE  OF  EACH  ELEMENT  OF  THE  FRONT  SIDE  OF  PLATE  (R) 
C  C  T( J, I, l J-TEMPERATURE  OF  EACH  ELEMENT  OF  THE  PLATE  BACK  SIDE  (R) 

C  C  HTCBAR( I,J)*CQNVECTIVE  HEAT  TRANSFER  COEFF.  AT  BACK  SURFACE 

C  OF  PLATE 

C  C  HTC(I»J)=HEAT  TRANSFER  COEFFICIENT  AT  FRONT  SURFACE  MATERIAL 

C  C  THT-ANGULAR  DISTANCE  FROM  ZERO-DEGREE  LINE  TO  ELEMENT  CENTER 

C  (OEG) 

C  C  R-RAO IAL  DISTANCE  FROM  PLATE  CENTER  TO  ELEMENT  CENTER  (FT) 

C  C  DELR-RADIAL  DIMENSION  OF  PLATE  ELEMENTS,  EQUALS  RA 

C  C  N-NUMBER  OF  RADIAL  ELEMENTS  USED 
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C  RAT»RADIUS  OF  JET  AT  THE  PLATE  (FT) 

C  VA-JET  mOCITY  AT  THE  PLATE  CF PS  > 

C  REFF-EFFECTIVE  RADIUS  OF  JET  (FTI 
C  OEFF-EFFECTIVE  DENSITY  OF  JET  (LB/CUOlC  FT) 

C  MEFF-JET  NACH  NO. 

C  N-LIHIT  ON  NUN8ER  OF  RADIAL  ELEMENTS  USED 

A  I-RADIAL  I  NOE  X 

A  J" ANGUL AR  INDEX 

A  KSTL "THERMAL  CONDUCTIVITY  OF  STEEL  (BTU/SEC-FT-F ) 

A  OSTL-DENSITY  OF  STEEL  (LB/CUBIC  FT) 

A  HC-SPECIFIC  HEAT  OF  STEEL  CBTU/LB-F) 

A  DA-OENSITY  OF  COMBUSTION  PRODUCTS  AT  2500F  (LB/CUBIC  FT) 

A  NJA-V1SC3SITY  OF  COMBUSTION  PRODUCTS  AT  2500F  (LB/FT-SEC) 

A  MBAR-THERMAL  CONC’IC  TI VITY  OF  COMBUSTION  PRODUCTS  AT  1500F 

(BTU/SEC-FT-F) 

A  MUBAR-VISCQSITY  OF  COMBUSTION  PRODUCTS  AT  1500F  (LB/FT-SEC) 
A  CPBAR-CONSTANT  PRES.  SPECIFIC  HEAT  OF  COMB.  PRODUCTS  AT 
1500F  (8TJ/L3-F) 

A  M-NUMBER  OF  ANGULAR  ELEMENTS  USED 

A  DELTHT" ANGUL AR  DIMENSION  OF  PLATE  ELEMENTS 
A  DAIR-DENSTIY  Oc  AIR  (LB/CUBIC  FT) 

A  PEFF-EFFECTIVE  PRESSURE  OF  JFT  (PSIA) 

A  RP-GAS  CONSTANT  FOR  PROPANE  I FT-L B/ LB M-R ) 

A  GC-DIMENSIONAL  CONSTANT  (L3M-SEC  SQ/LB-FT) 

A  GAM-SPECIFIC  HEATS  RATIO  OF  PROPANE 

A  SIGMA-STEFAN-BOLTZNANN  CONSTANT  ( BTU/ SEC-SQ FT • R**4 ) 
TEMP-ENTROPY  DIAGRAM  FOR  PROPANE.  P  -  ATM 
TK  «  DEG  K 

SF,SG  «  CAL/( MOLE-K) 

HF.HG  ■  CAL/MOLE 

VF#VG  «  CM3/M0LE 

PR  •  RESERVOIR  P RE SSUR E» P S I A 

XR  -  RESERVOIR  QUALITY*  0 IME NS  I ONL E S S 

MW  -  MOLECULAR  WEIGHT  OF  PROPANE*  GM/MOLE 


TEFFK 

MW  PC 

VEFFC 

ROEFFC 

AEFFC 

RAEFFC 

HEFFC 

UCPC 

HE  JC 

TCPK 

ROCPC 

HCPC 


TEMP  OF  PROPANE  AT  P-1ATM,  (K) 

MOLECULAR  WEIGHT  OF  PROPANE#  (GM/MOLE) 

VELOCITY  OF  EXPANOED  PROPANE  AT  1  ATM*  (CM/SEC) 

DENSITY  OF  EXPANDED  PROPANE  AT  1  ATM*  (GM/CMP*3) 

AREA  OF  EXPANDED  PROPANE  AT  1  ATM*  (CM**2> 

RADIUS  OF  EX0 ANDED  PROPANE  AT  1  ATM*  (CM) 

ENTHALPY  OF  EXPANDED  PROPCANE  AT  1  ATM#  (ERG/GM) 
INITIAL  VELOCITY  OF  COMB.  PROD.*  (CM/SEC) 

TOTAL  INTERNAL  ENERGY  OF  EXPANDED  PROPANE  DEPOSITED  TO 
COMBUSTION  PRODUCTS*  (ERGS/GM) 

INITIAL  TEMPERATURE  OF  COMB.  PROD.*  (K) 

INITIAL  DENSITY  OF  COMB.  PROD.*  (GM/CM**3> 

ENTHALPY  QF  COMB.  PROD.*  (ERGS/GM) 
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CPCPC 


NWCPC 

AOC  PC 

MACHCP 

ROAIRC 

ACPC 

RACPC 

AC 

VAC 

RAO  AC 

HEQCPC 

HAC 

TAK 


CONSTANT  PRESSURE  SPECIFIC  HEAT  OF  COMB.  PROD.# 

( ERG/( GM*K) > 

MOLECULAR  WEIGHT  OF  COMBUSTION  PRODUCTS#  (GM/MOLE) 
SPEED  OF  SOUNO  OF  COHB.  PROD.#  (CM/SEC) 

INITIAL  MACH  NO  OF  COHB.  PROD.#  (-) 

DENSITY  OF  AIR#  (GM/CM**3I 

INITIAL  AREA  OF  COMBUSTION  PRODUCTS#  (CM**2) 
INITIAL  RADIUS  OF  COMB.  PROD.#  (CM) 

DISTANCE  FROM  INITIAL  COMB.  REGION  TO  PLATE#  (CM) 
VELOCITY  OF  MIX  AT  PLATE#  (CM/SEC) 

RADIUS  OF  MIX  AT  PLATE#  (C") 

TOTAL  INTERNAL  ENERGY  CF  MIX  AT  PLATE#  (ERG/GM) 
ENTHALPY  OF  MIX  AT  PLATE 
TEMPERATURE  OF  MIX  AT  PLATE#  (K ) 


READ  ALL  INPUT  DATA 

1  READ(1» 1G00#END*2000) 

READ  (1#  6)  NNLOOP 
DO  8  1-1 » NNLOOP 
8  READ  (1# 1010)PPR(I)» M3F( I) 

5  READ(1#101C)XR*MW»RP 
DO  10  i»i,ie 

10  READ(1.1010)P(I)#TK(I)#SF(I)#SG(I)»HF(I)#HG(I)#VF(I)*VG(I) 
READ  FLAMf  DIMENSIONS 
READ  (1. 920)CD»RAT#TA 
READ  PLATE  DIMENSIONS 

READ  ( 1#  930 )  PLOST#PBDST#IGDST#BETA,THSTL#PLATE 
READ  INSULATION  DIMENSIONS 
READ  (1#910)THI,KI 
READ  SOOT  DIMENSIONS 
READ  ( l#  91 0 ) TH 3 »  KS 

READ  PROPANE  COMPRESSIBILITY#  LIMIT  OF  ITERATIONS 
READ  (1# 94 0)Z, LIMIT 
READ  TIME  INCREMENT 
READ  (1#945)DELT 

RE  AO  EMISSIVITY  OF  SURFACE  MATERIAL 
READ  ( 1 #  94 5 ) E 

READ  (1#  53) (THSAD( II )# II*1»125) 

DO  9997  II* l# NNLOOP 
DO  9997  JJ*1» NNLOOP 


DEFINE  CONSTANTS  FOR  PROPANE  PROBLEM 

GAM-1. IS 
RCONC-a. 31E7 
MRATIO-7 .6928 
HAIRC-5.023E10 


o  o  o 
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TAIRK-298.15 
SC-32.17 
DAIR-0.0765 
OELTHT-O. 5235987757 
MUA-3.5E-5 
DA-0.013 
KSTL-6.944E-3 
DS  TL -490 . 

HC-0.11 

CPBAR-0. 313 

MU3AR-2.76E-5 

KBAR-1.16E-5 

SIGNA-4.76111E-13 

SIGNAC-5 .67E-5 

EPSLNC-0.2 

PI-3.141592654 

UNAXC-1600.*12.*2.54 

RT«RAT*CD+*.5 

CALCULATE  RESERVOIR  CONDITIONS 

PR-PPRt II) / 14 • 7 
NRATIO*  HOF ( J J ) +2  . 

1-1 

15  IF(P(I)-PP)25»30.20 

20  I-I+l 
GO  TO  15 

25  ISAVE-I 

PR AT  1 0- ( PR-P( 1-1 ) ) /( P( I )  —  P  ( I-l  ) ) 

TR  -TK( I-l )+(TK(I»-TK< I-l) )*PRATIO 
SFR-SF(  1-1 )+<SF< 1 1 -SF ( I-1))*PRATI0 
SGR-SG(I-1)+(SG(I)-~G(I-1))*PRATI0 
HFR-HF( I-l )*<HF ( I ) -H F ( I-l) )*PRATIO 
HGR-HGt  I-l )*(HG( I)-HG< I-1))+PRATI0 
VFR-VFl  I-l  )  +  (  VF  ( I)-VP(  I-l)  )*PRATIO 
VGR«VG( I-l )♦( VG< I)-VG( I-1))*PRATI0 
SR  - ( l.-XR) *SFR+XR4SGR 
HR- < l.-XR )*HFR*XR*HGR 
VR«(1.-XR)*VFR+XR*VGR 
GO  TO  35 

30  ISAVE-I 
SFR-SF(  I ) 

SGR «$G(  I  ) 

HFR*HF(  I  ) 

VFR* VF( I  ) 

VGR- VG(  I  ) 

SR-(1.-XR)*SFR+XR*SGR 
HR  -  ( l.-XR)  ♦HFR*’XR*HGR 
VR«(1.-XR)*VFR+XR*VGR 
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35  WRITE (2, 1000) 

WRITE(2,1025)PR,XR,TR,SFR»SGR#SR#HFR#HGR,HR,VFR»VGR>VR 

ISENTROPIC  EX  PANS  I  ON  FROM  RESERVOIR  CONDITIONS  TO  CHOKED  FLOW 
AT  NOZZLE  THROAT 

WRITE (2# 1028) 

J-ISAVE 

K-l 

40  X( J)-CSR-SFCJ) )/(SG( JJ-SF ( J »  > 

VC  J ) " Cl •— X  C  J) ) *VFC  J)+X  C J ) *VG  C  J  ) 

H(J)»C1.-XCJ))*HFC J)*XCJ)*HGCJ> 

SCJ)-SR 

RHOC JI-MW/VCJ) 

VELCJ)*S0RT(2»*(HR“H(J))*4il86E7/HW) 

WDOTACJ)-RHQ(J)*VEL(J) 

WRITE C2,1030)XCJ), VC J > »H C J ), S C J  )»RHO C J )# VEL C J ) , WDOTA C J )* P C J 1 
IFCK.EQ.2)G0I0  45 
J«  J+l 
K-2 

GOTO  40 

45  IF  (WDOTA(J)-WDOTA(J-l)) 55,55, 50 
50  J-J+l 
GOTO  40 

CURVE  FIT  FOR  PT  AT  WDOTA  HAX 
FIT  PARABOLA  THROUGH  i,J-l  AND  J-2 

55  ENUNM  WDOTA  (J)-WODTA<  J-l ) )*( P C J  )-P C J-2 ) > -C WDOTAC J )-WDQTA C J-2 )) 
♦*CP ( J )-P ( J-l) ) 

DENOM«( P(J)**2-P< J-l )**2)*CPC J)-P( J-2) l- (P C J )**2-Pl J-2 )**2) 

♦*( P ( J )-P ( J-l) ) 

AA-ENUN/DENOM 

B»( (WDOTA(J)-WDOTA(J-l ) )-AA*  (P(J)*+2-P(J-l)**2))/CP(J)-PC  J-l ) ) 
PT«-B/(2.*AA) 

INTERPOLATION  FOR  THROAT  CONDITIONS  WHEN  WDOTA  IS  MAXIMUM 

WRITE  (2,1021)  R AT,C D, RT 
RATIO- (PT-PC  J-l))  / CP  CJ)-i»C  J-l )  ) 

TT  «TK(  J-l  H-CTKC  J)-TK(  J-l))* RAT  10 
SFT-SF(J-1)+(SF(J)-SF(J-1))+RATI0 
SGT»SG(J-1)MSG(J)-SS(J-1)>*RATI0 
HFT«HF<  J-l  )  +  (4F<  J)-HF(  J-1H  +  RATI0 
HGT«HG( J  — 1 ) ♦( HG  < J)-HS( J-l))* RAT  10 
VFT«VF< J-l )  +  (VF  C JJ-VFC J-l > )*RATIO 
VST-VGt J-l )+CVG( J)-VG( J-1))*RAT10 
XT*(SR-$FT)/(SGT-SFT) 

VT«C1.-XT)*VFT+XT+VGT 
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HT»(1.-XT)*HFT+XT*HGT 

RHQT-HN/VT 

VELT«S0RT(2.*<HR-HT)*4.186E7/MW> 

MDOTAM«RHOT*VELT 

NOOTH»2DOTAM*PI*<RT*12.*2.54)**2 

MRITEC2»1035)PT#XT,TT#SFT,SGT#  SR#HFT»HGT>HT, VFT, VGT*VT* RHOT*VELT# 

♦  MOOT AH» K  DOT H 

ISENTROPIC  EXPANSION  OF  PROPANE  FROM  NOZZLE  TO  ATMOSPHERIC  PRESSURE 

PEFF-1.0 

1-1 

60  IF ( P ( I> •EO. PE FF) GO TO  65 

i»i*i 

GOTO  63 
65  1SAVE-I 

TEFPK«TK(I  ) 

XEFF-<SR-SF<ISAV6) >/(SG< ISAVE )-SF(ISAVE» ) 

S  P^OL  ■(  1  •  - XE F  F  I  •  Vr  ( I  S  A  V£  )  ♦XE F  F ♦  V6(  I S  AVE) 
ENTH“(I«“XFFF)*HF(ISAVE)*XEFF*HGIISAVE) 

DEMS-fN/SPVOL 

FI JX-0ENS*S3«T( 2.*<HR-ENTHI*6.186E7/MW> 

AEFF.MOOTN/FL'JX 

RFFF.7)RT( AEFF/ PI  1/(2* 56*12.  ) 

XeFF-(H.UX/DcNS)/(2.5<i*12.) 

TE^F-TA 

NUNtRICAL  OIFFtRENTIATION  OF  T-S  DIAGRAM  FOR  PROPANE  AT  S«SR  AND  PEFF 
TO  OBTAIN  SPFED  OF  SOUND 

Pl.PEFF^O. 1 
I-l 

70  IFIPm.LT.PDGOTO  75 
I  -  t  ♦  I 
GJTO  70 

75  RATI0«(PI-P(I-1 )  )/ (p (I  l-P(I-l)  ) 

SF1«SF( I-1)*(SF( II-SFC I-1»)*RATI0 
SG1«SG(  1-1  >  ♦  (  S  G  C  I)-SGU-l»  »*RATIO 
XI«(SR-SFU/(SG1-SF1) 

VFl-VFI I-l )♦< tfF ( I)-VF( 1-1) )*RATIO 
VGl-VGI  I-l  )♦(  VG(I)-VG(  I-m*RATIO 
VI«(1.-X1)*VFI*X1*VG1 
RHOl-NW/Vl 
P2«PEFF-0.1 
1-1 

80  IF(P(I).LT.P2)SUT0  35 
I-I  +  l 

goto  eo 

85  RATIO«(  P2-P(I-mMP(  I)-P<  I-i  >  ) 
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SF2-SF<I-l)*<SFm-$F<I-l)l*RATI0 
$G2-SG(I-1)*<SG(I)-SG<I-1>>*RATI0 
X2« ( SR-SF2  )  /( SG2-SF2 ) 

VF2»VFCI-l)+CVFCII-VF(I-l)i*RATI0 

VG2-VG<I-ll*(VGm-VG<I-l))*RATI0 

V2-(1.-X2)*VF2*X2*VG2 

RHQ2-MW/V2 

DRHODP-ABS CRH02-R401)/ (2. *0.1) 

S0NIC«S9RT<  <14.7*4  54.*930./(2.54**2>  )/DRHQDP>/ C12.42. 54) 
NEFF-VEFF/SONIC 

RH0EFF-144.*14.7*PEFFMZ*RP*TEFF> 

PEFF«14.7 

PEFFC«PEFF*6.895E4 

0EFF-RH0EFF 

WRITE (2* 1040> REFF. VEFF»TEFFf SONIC* KEFF#RHQEFF# PE FF 
MWPC-MW 

VEFFC-VEFF*2. 54*12. 

R0EFFC-PEFFC*MWPC/(RC3NC*TEFFK) 

AEFFC-AEFF 

RAEFFC»REFF*2. 54*12 

HEFFC-ENTH*(1/MWPC)*4.186E7 

WRITE (2# 1045)PEFFC*VEFFC*  ROE  FFC#AEFFC»RAEFFC#HEFFC#TEFFK 
EMO  OF  EXPANSION  TO  ATM  CONOITIONS 


BEGIN  COMBUSTION 


UMIXC«( l.+MRATIO) ) *VEFFC 
C  HEQC-HEFFC*UVEFFC**2.-UCPC**2.)/2.) 
TU0X0K«133.9*MRATI0*337.3 
DTDXMK*-2»  4734*MRATIO*50.73 
0TDUMK--0.6878+MRATI 0*14.091 
TUOXXK*TUOXOK+OTOXMK*XEFF 
TCPK»<JMIXC/UMAXC>*DTDUMK*TUOXXK 
MWC PC »0. 74 5 6* MR  ATI  0*17 .44 
HN0T-34786.0/MWCPC 
DH0XM«-0.99**1RATI3*19.42 
HX0M*7.Q*M RAT  1 0-129.7 
HCP=HXOM+DHDXM*XEFF 
HCPC»MCP+HN0T)*4. 186E7 
ROCPC»( PEFFC*MWCPC )/ (RC0NC4TCPK) 
GAMCPC=-8.254E-3*MRATT0* 1.3635 
UCPC=UMI XC 

AOCPC«(GAMCPC*RCONC*TCPK/MWCPC  >**0.5 

CPCPC*1 . 53  BE7 

MACHCP*UC°C/AQCPC 

WRITE(2»1055)TU0X3K»0TDXMK» 0  TDUMK»  TUOXXK*TCPK 

WRITE(2»1056)NN0T»040XM,HX0M,HCP#HCPC 

WRITE(2»1057)R0CPC»GAMC9C*UCPC#A0CPC#CPCPC 
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ROAIRC«0AIR*1.6O5E-2 

ACPC»CROEFFC*AEFFC*VEFFC )*<l.*MRATIO)/ (RQCPC*UCPC) 
RACPC-SQRTUCPC/PI) 

WRITE(2»1065)ROAIRC#ACPC#  RAC  PC 
A« ( PBOST-I GOST ) 

AC«A*12.*2.54 

IFCMACHCP.LT. l.OlSOTO  90 

KAPPA»0.063*(MACHCP**2.-1)**(-0»15) 

SOTO  95 

90  KAPPA «0. 08 *{1.-0»16*HACHCP)*(R0AIRC/R0CPC)** (-0.22) 

95  DENOM«(KAPPA*CAC/RACPC )* (ROAIRC /ROCPC ) **0. 5 )-0. 7 
OEMOMl-C  C(0EN0M+.7)/KAPPA)*.102)-.7 
WRITE (2» 1060) AC »KAPPA»DEN0M 
VAC»UCPC*( l.-EXP (-1. /DENOM) ) 

VA-VAC/C2. 54*12) 

RADAC«CAC*0.07)-0.70*RACPC*( ROAIRC /ROC  PC >**(-0.50) 

RA-RADAC/C2. 54*12) 

HEQCPC«HO>C-MUCPC**2.  )  /2. -HA IRC 

HAC»HAIRC-(VAC**2. )/ 2. *HE QC PC *<1. -EX  PC-1 /DEN0M1) ) 

TAK-TCPK  +  CHAC-HCPC  )/CPCpC 

PO-PR 

T0«  TR 

00“ MW/ VR 

WRITEC2#62)P0»T0#00»RT#TA»A#BETA#THSTL# PLATE 
WRITE(2»63)THI>KI,THS»KS»Z»LIMIT#E 

WRITEC2#965)PEFF#TEFF,MEFF#DEFF#AEFF#REFF#VEFF#VA#RA 
WRITEC2»1067)V4C»RA0AC»HEQCPC#HAC#TAK 
Z00TC-(PLDST-P80ST)*2. 54*12.0 
UZWPC-/AC*(l-EXP(-2*(Z0'JTC/RACpC)**2.) ) 

WRITEC2*  1022)  ZOdTCuJZHPC 

VASTPCC II# JJ)«UZWPC*C l.M ROCPC /ROAIRC-1. )*UZWP C/UCPC )*♦. 5 
DELR-RA 

N-PLAT=/0ELR*0.5 
IF  ( N.GT .3 ) GOTO  301 
N-4 

OELR-PLATE/CH-O.S) 

301  WRITEC2#52)0ELR#N 

CALCULATE  HEAT  TRANSFER  COEFICIENTS  AT  FRONT  SURFACE  OF  PLATE 
M«7 

205  BETA«BETA*PI/180.0 
RCEN«.2*DELR 

HTC(1»1)*C0.04Q*KBAR/RCEN)*CCP8AR* MU BAR/KB AR)**0.4 
♦*C0A*VA*RCEN/MUA)**0.9 
HTCC-HTC(1#1)*2.044£7 

300  DEG  K  IS  THE  ASSUMED  WATER  TEMPERATURE  OF  WATER  IN  THE  PROBE 
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QCIII*  JJ)»HTCC4(TAK-300. )+SIGMAC4EPSLNC4CTAK*44-300.*44.  ) 

1-1 

WRITE (2*  49 1 
WRITE (2* 48  ) 

DO  206  J  *1#  7 

206  HTC(1#J)-HTC(1,1) 

WRITE(2#950)(HTC(I*J )»  J»1>N), I 
DO  22C  I  *2»N 
RCEN-(I-1I*DELR 

DO  210  J  *1 » M 
TNTCEN»DELTHT*(J42) 

IP(J-4)207»208,209 

207  HTC(I*J ) *0*0404 ( KB  AR/RCEN l*( C P BAR* HU BAR/ KB AR 1*40 • 44 
♦((DA4VA4RCEN/MUA)4(l.O*C0S(BETA) >4 < 1 +C 0$ ( BET  A) 4$ IN < THTCE N  )  > ) 4*0. 9 

GOTO  210 

208  HTC(I*J)-0.0404(KBAR/RCEN)4(CPBAR4HUBAR/KBAR>440.44 

♦  ( ( 0A4VA4RC  EN/NUA ) 1 440*  9 
GOTO  210 

209  HTC (I»J) «0.0404(K3AR/RCEN)4( CPB AR4MUBAR/ KB AR >440.44 

+  ( <DA4VA*RCEN/*1UA  )*(1. 0-COS  (BETA  )  )4<1+C0S  (BET A>4$ IN(THTCEN) > )4*0. 9 

210  CONTINUE 

WRITE(2»950)(HTC(I#J  )* J»1,M)#I 
220  CONTINUE 
WR ITE (2*  7 ) 

9997  CONTINUE 

DO  9996  I* 1* NNL OOP 
DO  9996  J* l.NNLOOP 
WRITE (4» 3001)VASTPC(I# J) 

9996  WRITE(4»3001)QC(I#J> 

GOTO  2000 

6  FORMAT  (110) 

3001  FORMAT  (E10.4) 

7  FORMAT  ( 1H 1 5 
910  FORMAT  (2F20.5) 

920  FORMAT  (3F20.5) 

930  FORMAT  (6F10.5) 

940  FORMAT  (F20.5.I3) 

945  FORMAT  (F10.4) 

48  FORMAT  (8X*'9G  DEG  '*8X»*120  DEG  *»7X»'150  DEG  '»7X»»180  DEG  • #7X* 

♦  •210  DEG  •  1 7X»  *240  OEG  '#7X#*27C  DEG  '#8X,»R»,//> 

49  FORMAT  (H1»40X,MTC  MATRIX  (  BTU/SEC-SOFT-F)  •* // ) 

950  FORMAT  ( 7F 1 5. 7* 5 X . 1 3# •  OELR»»/> 

52  FORMAT  CODEIR  «  • »  F14. 6#  1 2X*  •  N  »'#I4 ,/) 

53  FORMAT  ( 80A1> /# 45A 1 J 

62  FORMAT  ('IPO  «• * E16.4# 10X# 'TO  ■  »#E16.4#10X#  •  DO  ■*»E16*4#  10X»  '  R  J  •' 

♦  »  E 16. 4» / /»  •  TA«'»E16.4#10X»*A  «'»E16*4*10X»,BETA  ••,E16.4#10X, 'THS 

♦  TL  •  *  »E 16*  4#//# •  PLATE  «'#E16.4W/> 

63  FORMAT  (*OTHI  « • ,E16.4,10X, • KI  •*» E16.4» 1CX# *THS  » E16. 4#10X, *KS 

♦  «'»E16.4,//f  •  Z»'*E16.4» 10X, *IIMIT  • I10.10X, •£  «',F5.3,//) 


965  FORMAT  ( 6H  PEFF-#E16.4#/» 6H  TE FF»# E 16. 4# /# 6H  MEF F»#16. 4# /#6H  OEFF- 
♦#E16.4#/#6H  AEFF-»E16.4#/#6H  REFF-# E 16 . 4, /, 6H  VEF F-#E16.4#/,4H  VA- 
♦#E16. 4# / #4H  R A-  #  E 1 6.  4#  //  ) 

87  FORMAT  ( *ODELT  - • > El  2. 4# 10X# • DELT1  « •# El  2.4# 10X# • DELT2  -'#£12.4,10 

♦  X# • TRES  *  *  #  E12. 4 ) 

350  FORMAT  (•©♦♦♦♦^  TIME  INCREMENT  CHOSEN  IS  TOO  LARGE#  MUST  BE  LESS  T 
♦HAN  OR  EQUAL  T0»#E16.4#*  SECONDS^****' ) 

487  FORMAT  ( 1H1#125A1# /) 

489  FORMAT  (•  •) 

850  FORMAT ( 1H  > F14. 2# F14. 3# 8X» • C ENTER' # E 20 .4# F 16 .2# F17.2 ) 

900  FORMAT  (1H  #F14. 2# FI 4. 3# F14. 1# E20. 4# F16. 2# F17. 2) 

985  FORMAT  ( IX# 99( •♦•)#//) 

1000  FORMAT  (80H1 

♦  > 

1010  FORMAT! 6E12.6) 

1025  F3RMAT<1H0»11X»20HRESERV0IR  CONDITIONS/ 

♦  2X#7HPR(ATM)#22X#E9.3/ 

♦  2X#11HXR(QUALITY)#18X#E10.4/ 

♦  2X#5HTR(K)#24X#E11.5/ 

♦  2X#15HSFR(CAL/M0LE-K)#14X#E1C.4/ 

♦  2X#15HSGR(CAL/M0L£**K)#14X#E1C.4/ 

♦  2X,15HSR.(CAL/N0LE-K)#14X#EI0.4/ 

♦  2X#13HHFR(CAL/M0LE)#16X#E10.4/ 

♦  2X#13HHGR(CAL/NOLE)#16X#E10.4/ 

♦  2X#13HHR  (CAL/M0LE)#16X#E10.4/ 

♦  2X#13HVFR( CM3/10LE ) ,16X#E10. 4/ 

♦  2X#13HVGR(CM3/inLE)#16X#E10*4/ 

♦  2X#1 3HVR  (CM3/M0LE),16X#E10.4I 

1U21  FORMAT  C 1H0#2X# 23HRADIUS  OF  THROAT  ( FEET )# 9X#E10 .4/ 

♦  2X.6HCD  (-  )#  26X#  E10.4/ 

♦  2X#  30HEFFECT IVE  THROAT  RADIUS  (FEET) #2X#E10.4> 

1028  FORMAT (  HO#  10X,  32 HSE  ARCH  FOR  CHOKED  FLOW  AT  THROAT/ 

♦  5X#1HX#5X,8HV(CM3/M)#2X#8HH(CAL/M)#2X»  8HS( C/ M-K )» 2X» 9HRH0( G/CC )# 

♦  1X#8HVEL(C/S)#4X#4HFLUX#6X#6HP(ATM>) 

1030  FORMAT ( FB. 4#3X# F8. 1# 3X#F 7. 1# 3X»F6. 2#  3X# F 8 . 5# 2X# F 8 . 1# 2X» F7.2# 4X # 

♦  F5.1 ) 

1035  FORMAT(1HO#11X#17HTHROAT  CONDITIONS/ 

♦  2X#7HPT(ATM)#22X»E9.3/ 

♦  2X#11HXT(QUALITT)#18X#E10*4/ 

♦  2X#5HTT(K)#24X#E11.5/ 

♦  2X#15HSFT(CAL/M0LE-K)#14X#E1G.4/ 

♦  2X# 15HSGT( CAL /M OL E-K) #14X»E 10.4/ 

♦  2X#l 5HS T  (CAL/MQLE-K)#14X»E10.4/ 

♦  2X#13HHFT(CAL/M0LE)#16X#E10.4/ 

♦  2X#13HHGT(CAL/M0LE)#16X#E10.4/ 

♦  2X#1 3HHT  (CAL/M0LE)#16X#E10.4/ 

♦  2X#13HVFT(CM3/M0LE)#16X#E1C.4/ 

♦  2X#13HVGT(CM3/M0LE)#16X#£10.4/ 

♦  2X#13HVT  (CM3/M0LE)#15X#E10.4/ 


CO 


4  2X»12HRH0T(GM/CM3 )>17X»E10«4/ 

♦  2X*L2HVELT (CM/SEC ) * 17X#E10» 4/ 

4  2X#17HFLUXT(GN/CM2-SEC)#2X#E10.4/ 

4  2X*13HFL0WT (GM/SEC)*4X»E10*4) 

1040  FQRMAT(1H0#7X»20HEFFECTIVE  CONDITIONS/ 

4  2X*8HREFF( FT  >  » 1 2X  »  E 10 . 4/ 

4  2X»9HVEFF(FRS)»11X»E10.4/ 

♦  2X#7HTEFF(R)#13X#E10»4/ 

4  2X»10HS0NIC ( F PS ) » 10X* E10*  4/ 

4  2X#4HMEFF#16X»E10.4/ 

4  2X*14HRH0EFF(L8/FT3)#6X#E10.4/ 

4  2X»10HPEFF(PSIA)»10X>E10.4) 

1045  FORMAT (1H0»7X»31HEFFECTIVE  CONDITIONSCCGS  UNITS)#/* 

4  2X»5HPEFFC»15X»E10.4/ 

4  2X#5HVEFFC*15X»E10.4/ 

4  2X,6HR0EFFC»14X,E10.4/ 

4  2X»5HA£FFC#15X/E10. 4/ 

4  2X#6HRAEFFC,14X»E10.4/ 

4  2X#5HHEFFC#15X»E10.4/ 

4  2X,5HTEFFK#15X#E10.4) 

1055  FORMATC 1H0# /# 2X# lOHTUOXOK  <K ) , 6X#E 10.4/ 

4  2X# 10HDTDXMK  (K >* 6X*E 10.4/ 

4  2X» 13HDTDUMK  (K ), 6X,E 10.4/ 

4  2X# 10HTUJXXK  (K )#6X#E10«4/ 

4  2X#  8-1TC PK  <K)#8X#E1G.4> 

1056  F0RMAT(1H0»/#2X#13HHN0T  ( CAL /GN ) #3X, E10. 4/ 

4  2X» 144DH0XM  <CAL/GM)#2X#E10.4/ 

4  2X» 13HHX0M  (CAL/GM)#3X#E10.4/ 

4  2X » 12HHCP  ( C AL/GM  )  #4X*E10*4 / 

4  2X#13HHCPC  (ERG/GM)#3X#E10. 4) 

1057  FORMAT( 1H0» /, 2X, 13HR0CPC  (GM/CC  )#3 X* E1C. 4/ 

4  2X.10HGAMCPC  (- ), 6X,E10.4/ 

4  2X» 13HUCPC  (CM/SEC  )*3X#  E10*  4/ 

4  2X» 14HAQC  PC  (CM/SEC)#  2X*E10.4/ 

4  2X» 16HCPCPC  (ERG/GM  K),1X,E10.4> 

1022  F0RMAT(1H0»2X#334DIST4NCE  FROM  PLATE  TO  PROBE  (CM  )# 2X» E10.4/ 
4  2X»  31HVEL  AT  PROBE  WITH  PL ATE (C M/ SEC )# 4X, E10.4) 

1060  FORMAT( 1H0» /»  2X » 24 AC » 1 8X# E10. 4 / 

4  2<# 5HKAPPA#15X,E10.4/ 

4  2X#5HDEN0M#15X#clC.4) 

1065  FORMAT! 1H0, /, 2X# 6HR0AIRC#14X, E10.4/ 

4  2<»4HACPC»16X»E19.4/ 

4  2X*54RACPC*15X,10.4) 

1067  FORMATl 1H  * /, 2X, 34/AC# 17X, E1G. 4 / 

4  2X#54RADAC»15X,E10.4/ 

4  2X»6HHEQ:PC#14X#E1C.4/ 

4  2X#34HAC»17X, E1C.4/ 

4  2X»34TAK*17X#E10.4> 

2000  CALL  EXIT 
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Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet ,  fold  as  indicated,  staple  or  tape  closed,  and  place 
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2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
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